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Abstract 

A  statistically-based  forecasting  tool  is  developed  for  Dover  AFB,  McGuire  AFB, 
and  Andrews  AFB  for  dissipation  times  of  fog  and  low  stratus.  Probability  forecasts  are 
produced  at  hourly  increments  from  0-6  hours  for  the  most  extreme  reductions  in 
visibility  (less  than  0.5  mi)  and  ceilings  (below  200  ft).  Forecasts  are  based  on  surface 
observations,  upper  air  observations,  and  climatological  parameters. 

Ceiling  forecasts  at  Dover  AFB  and  McGuire  AFB  show  improvements  over 
conditional  climatology  ranging  from  1-51%  with  an  average  improvement  of  19.2% 
when  verified  against  an  independent  data  set.  McGuire  AFB  visibility  forecasts  show  an 
average  improvement  over  conditional  climatology  of  3%.  These  findings  are  of 
particular  importance  to  the  Air  Force  in  general  and  specifically  to  the  15**^  Operational 
Weather  Squadron  (15  OWS)  who  produces  forecasts  for  these  airfields.  Implementing  a 
method  superior  to  conditional  climatology  is  expected  to  provide  improved  forecasts  and 
flight  operations  for  these  sites. 

The  two  forecasts  for  Andrews  AFB  show  relatively  low  mean  square  errors,  but 
are  unable  to  consistently  improve  on  conditional  climatology,  demonstrating  an  average 
decrease  in  forecasting  skill  of  42%.  Small  samples  of  data  could  be  the  reason  for  the 
decrease  in  skill.  The  Dover  visibility  forecast  also  shows  negative  forecast  skill,  with  an 
average  decrease  of  39%. 

The  method  is  a  success  in  producing  forecasts  for  ceiling  and  visibility  criteria 
that  had  never  previously  been  examined.  Further  research  on  the  technique  could 
produce  a  powerful  tool  consistently  able  to  defeat  conditional  climatology.  It  is 
suggested  that  the  15  OWS  incorporate  this  methodology  into  their  operational 
forecasting  routine. 
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A  STATISTICALLY  BASED  METHOD  EOR  PREDICTING  EOG  AND  STRATUS 


DISSIPATION 

I.  Introduction 

Weather  foreeasters  in  general  are  unable  to  consistently  beat  persistence  in  short¬ 
term  (0-6  hour)  forecasts  of  ceiling  and  visibility.  Dagostraro  et  al.  (1995)  and  Dallavalle 
and  Dagostraro  (1995),  showed  that  not  only  is  persistence  competitive  with  National 
Weather  Service  (NWS)  forecasts,  but  that  NWS  forecasts  have  not  improved  in  the 
period  from  1985-1995.  The  same  inadequacies  found  in  the  NWS  are  experienced 
throughout  Air  Eorce  (AE)  weather.  The  standard  to  beat  in  the  AE  is  conditional 
climatology,  because  it  enables  forecasters  to  incorporate  knowledge  of  what  has 
happened  in  similar  situations  in  the  past,  thus  it  is  considered  a  better  predictor  of  future 
weather  than  persistence  (Hilliker  and  Eritsch  1999). 

A  statistical  technique  based  on  Vislocky  and  Eritsch  (1997)  and  Hilliker  and 
Eritsch’ s  (1999)  work  could  provide  a  tool  for  the  15*’^  Operational  Weather  Squadron  (15 
OWS)  to  improve  forecasts  from  conditional  climatology  tables.  This  research  develops 
an  automated  statistical  data  system  using  an  observations-based  network  for  operational 
implementation. 

Conditional  climatology  combines  persistence  and  climatology  by  incorporating 
knowledge  of  outcomes  of  similar  weather  situations  in  the  past  in  order  to  provide  short- 


1 


term  forecasts  of  current  weather  parameters.  Because  it  represents  averaged  conditions 
over  a  long  period  of  time,  it  is  a  more  difficult  tool  to  beat  than  persistence  (Vislocky 
and  Fritsch  1997).  However,  conditional  climatology  tables  do  exhibit  some  problems. 
Murphy  and  Katz  (1995)  illustrate  that  small  sample  sets  in  certain  cells  can  lead  to 
unstable  probabilities,  as  would  be  especially  true  for  the  extreme  cases  examined  in  this 
research.  It  is  therefore  imperative  that  the  AF  develops  a  more  skilled  forecasting 
method  for  fog  and  stratus  dissipation. 

1.1  Statement  of  the  Problem 

Accurate  predictions  of  fog  and  low  stratus  dissipation  (for  the  purpose  of  this 
paper  low  stratus  is  defined  as  clouds  at  or  below  200  ft  AGL  and  the  term  fog  refers  to 
fog  and/or  low  stratus)  are  essential  to  all  aircraft  operations.  The  15  OWS  is  responsible 
for  producing  forecasts  for  1 1  airfields  throughout  the  northeastern  quadrant  of  the 
United  States,  stretching  from  North  Dakota  to  Washington  DC  (Fig.  1).  Aircraft 
supported  by  the  15  OWS  range  from  large  fixed  wing  aircraft  to  rotary  wing.  Among 
these  are  transport  aircraft,  which  have  critical  thresholds  for  flight  cancellations  of  200  ft 
ceilings  and  0.5  mile  visibility.  Accurate  ceiling  and  visibility  forecasts  are  critical  to 
ensure  personnel  and  essential  cargo  is  transferred  in  a  timely  manner.  Every  hour  of 
inaccuracy  costs  the  United  States  government  thousands  of  dollars  and  possibly  lives 
(FAS  2003). 

There  are  numerous  reasons  for  the  lack  of  skill  in  predicting  fog  dissipation.  Fog 
dissipation  is  highly  variable  situation  and  based  on  small  changes  in  a  variety  of 
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FIG.  1.  Location  and  ICAOs  of  the  15  OWS  Airfields.  Airfields  with  the  (♦)  symbol 
are  the  eritieal  bases,  home  to  aireraft  with  minimums  of  200ft/0.5  mi.  Adapted  from  the 
15  OWS  homepage  (15  OWS  2003). 


meteorologieal  parameters  (Anthis  &  Cracknell  1998).  Furthermore,  the  physical 
environment  presents  large,  localized  influences  sueh  as  snow  eover  (Johnson  1978), 
variable  terrain  and  bodies  of  water  that  eontribute  to  variability  in  dissipation  rates 
(Weiss  and  Gurka  1975).  Other  important  factors  that  would  aid  in  prediction,  such  as 
depth  of  the  fog  and  height  of  the  inversion  layer  are  not  routinely  measured  nor  easily 
obtained.  Vislocky  and  Fritseh  (1997)  noted  that  the  single  most  important  faetor  in 
short-term  foreeasting  is  the  latest  surfaee  observation.  For  this  reason  alone,  numerieal 
models  are  at  a  disadvantage  when  it  eomes  to  predieting  dissipation. 

Of  these  foreeast  diffieulties,  the  vast  differenees  in  terrain,  elimatology,  and 
relationship  to  bodies  of  water  among  the  airfields  in  the  15  OWS  Area  of  Responsibility 
(AOR)  represent  a  significant  challenges.  These  variations  present  foreeasters  with  a 
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TABLE  1.  Annual  Number  of  Days  with  Reduced  Ceiling  and  Visibility.  The  three 
highlighted  bases  are  home  to  transport  aircraft  with  flight  minimums  of  200  ft/0.5  mi. 
(15  OWS  2003). 


Ceiling/Visibility  Threshold  Criteria  (ft/mi) 


Airfield 

ICAO 

vis  <  7 

<3000/3 

<1500/3 

<1000/2 

<500/1.5 

<200/.5 

Minot  AFB  (ND) 

KMIB 

71 

62 

40 

26 

15 

4 

Grand  Forks  AFB  (ND) 

KRDR 

72 

66 

40 

26 

15 

4 

Ellsworth  AFB  (SD) 

KRCA 

47 

47 

33 

29 

18 

7 

Offutt  AFB  (NE) 

KOFF 

109 

73 

47 

33 

15 

4 

Scott  AFB  (IE) 

KDLV 

197 

69 

40 

26 

11 

4 

Wright-Patterson  AFB  (OFl) 

KFFO 

193 

84 

47 

29 

11 

1 

Fort  Drum  (NY) 

KGTB 

117 

88 

47 

33 

18 

4 

McGuire  AFB  (NJ) 

KWRI 

192 

77 

58 

44 

26 

8 

Fort  Belvoir  (VA) 

KDAA 

170 

69 

55 

37 

22 

7 

Andrews  AFB  (VA) 

KADW 

170 

69 

55 

37 

22 

7 

Dover  AFB  (DE) 

KDOV 

187 

58 

44 

29 

15 

4 

variety  of  challenges  from  different  types  of  fog  and/or  physical  mechanisms  contributing 
to  formation  and  dissipation.  Table  1  quantifies  some  of  these  difficulties  by 
summarizing  the  number  of  days  in  which  fog  impacts  operations.  The  data  illustrate  that 
low  ceilings  and  visibility  are  common  occurrences  throughout  the  AOR.  The  eastern 
seaboard,  home  of  a  majority  of  the  cargo  aircraft  in  the  AF  inventory,  experiences 
reduced  ceilings  and  visibility  on  average  half  of  the  days  of  the  year,  significantly 
impacting  military  operations. 


1.2  Research  Objectives 


The  overall  objective  of  this  research  is  to  produce  an  automated  ceiling  and  visibility 
forecasting  tool  for  the  15  OWS  that  will  improve  on  all  currently  employed  methods. 
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The  specific  research  goals  are  as  follows: 

1 .  Develop  a  statistically-based  system  using  surface  and  upper  air  observations 
as  well  as  climatological  predictors  to  accurately  forecast  fog  dissipation. 

2.  Completely  automate  system  to  ensure  ease  of  use. 

3.  Present  a  viable  forecasting  tool  for  the  15**'  OWS. 

The  next  topic  addressed  is  a  summary  of  past  research  on  successful  prediction 
methods  of  fog  dissipation.  Various  techniques  are  examined  including  numerical 
weather  prediction,  conditional  climatology,  observations-based  methods,  satellite 
imagery,  and  multi-source  methods,  which  combine  two  or  more  of  the  above  to  produce 
a  dissipation  forecast.  The  third  chapter  focuses  on  the  methods  employed  in  this 
research  to  develop  the  forecasting  tool  for  the  15  OWS.  The  fourth  chapter  provides 
detailed  results  of  the  research,  and  comparisons  to  conditional  climatology  are 
examined.  Finally,  conclusions  are  made  which  demonstrate  the  effectiveness  of  the 
newly  developed  method,  and  suggestions  are  offered  for  possible  improvements  in  the 
future. 

Using  statistical  techniques  for  forecasting  weather  is  not  a  new  concept, 
however,  with  the  development  of  physically-based  numerical  models,  statistical 
techniques  have  become  less  prevalent  in  short-term  forecasting.  Chapter  two  examines 
the  success  rate  that  statistically-based  forecasting  has  produced  in  the  past  for  short-term 
prediction  of  fog  dissipation. 
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II.  Literature  Review 


This  research  is  focused  on  the  dissipation  of  fog  and  does  not  address  fog 
formation  or  prediction.  This  section  focuses  on  several  different  methods  to  produce  fog 
dissipation.  While  it  focuses  mainly  on  statistically-based  methods  to  predict  dissipation, 
other  techniques,  such  as  numerical  models,  conditional  climatology,  satellite  imagery, 
and  multi-source  methods  are  examined.  Advantages  and  disadvantages  are  examined,  as 
well  as  areas  where  these  techniques  can  be  improved. 

2.1  Numerical  Models 

Numerical  models  are  favored  by  forecasters  due  to  their  ease  of  use,  but  are  often 
not  the  best  available  asset  for  short-range  forecasts.  Statistically-based  products  derived 
from  high-resolution  mesoscale  models  can  predict  fog  dissipation  times,  but  have  not 
proven  reliable.  Porter  (1995)  illustrated  the  fact  that  persistence  is  superior  to  mesoscale 
numerical  model  predictions  in  the  short-term.  Furthermore,  Vislocky  and  Fritsch  (1997) 
showed  that  while  Model  Output  Statistics  (MOS)  can  consistently  beat  persistence  and 
conditional  climatology,  there  are  still  better  methods.  The  NWS  (1981,  1995)  concluded 
that  observations  are  the  most  important  ingredient  in  producing  short-range  forecasts, 
leading  Vislocky  and  Fritsch  (1997)  to  suggest  that  a  surface-based  observing  network 
would  be  more  effective  than  MOS  guidance. 
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2.2  Conditional  Climatology 


As  mentioned  earlier,  eonditional  elimatology  offers  an  advantage  over  persistenee 
in  that  it  ineorporates  knowledge  of  what  has  happened  in  past  situations  to  produee  an 
aeeurate  probability  forecast.  Data  are  interpreted  from  tables  which  are  stratified  by 
station,  month,  time  of  day,  and  wind  direction  and  offer  probabilities  out  to  48  hours  on 
the  occurrence  of  select  ceiling  and  visibility  criteria.  Tables  are  applied  based  on  the 
current  ceiling  criteria  already  being  met.  As  mentioned,  these  tables  are  stratified  by 
time  of  day  and  wind  direction,  therefore,  the  only  required  meteorological  parameter 
needed  to  produce  a  forecast  is  the  current  surface  direction  at  the  station. 

This  methodology  is  favored  by  the  Air  Force  for  its  simplicity  and  its  consistent 
ability  to  produce  fairly  accurate  forecasts.  The  general  problem  with  any  climatological 
model  is  that  it  does  not  take  into  account  the  current  meteorological  event  occurring  at 
and  around  the  station.  Observations-based  networks  include  this  data  in  statistical 
methodology  to  improve  on  conditional  climatology. 

2. 3  Observations-Based  Statistical  Methods 

Enger  et  al.  (1964)  introduced  the  concept  of  statistical  techniques  for  2-7  hr 
prediction  of  ceiling  and  visibility  in  the  1960’s.  Experiments  in  Chicago  and  major 
cities  on  the  west  coast  of  the  US  compared  objective  forecasts  to  conditional 
climatology  and  subjective  forecasts.  Enger  et  al.  (1964)  used  450  binary  predictors  to 
forecast  future  values  of  36  ceiling  and  visibility  predictands.  The  results  showed  that 
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statistical  techniques  were  superior  to  conditional  climatology,  persistence,  and 
subjective  forecasts.  Improvement  ranged  from  -0.4%  to  33.0%  with  average 
improvement  of  19.2%  for  eeiling  foreeasts  and  12.8%  for  visibility  forecasts  (Enger  et 
al.  1964). 

Vislocky  and  Fritsch  (1997)  showed  more  reeently  that  observations-based 
statistieal  systems  are  more  effective  than  MOS  and  eonditional  climatology  for  short¬ 
term  ceiling  and  visibility  foreeasts.  Their  research  focused  on  25  major  cities  along  the 
eastern  eorridor  of  the  US  and  used  surfaee  observations  from  a  predetermined  area 
around  each  airfield  to  foreeast  future  conditions  of  ceiling  and  visibility  (Visloeky  and 
Fritseh  1997).  Hilliker  and  Fritseh  (1999)  proved  that  a  similar  method  eould  be  used  to 
predict  fog  dissipation  at  the  San  Francisco  airport. 

Vislocky  and  Fritsch  (1997)  used  a  number  of  predietands  (Table  2)  based  on 
signifieant  airfield  flight  restrictions.  They  embarked  on  a  pilot  study  to  determine  the 
optimal  number  of  surrounding  stations  and  predictors  to  use  in  order  to  eliminate  data 
saturation.  The  results  produced  five  significant  findings: 

1 .  The  optimal  number  of  stations  to  eonsider  inereases  with  desired  lead-time 
and  is  not  constant.  Variations  exist  with  season,  latitude-longitude,  and 
weather  among  other  faetors. 

2.  The  ceiling  and  visibility  predietors  were  the  most  important. 

3.  Non-predietand  observations  outside  the  optimal  number  of  stations  had  very 
little  impaet  on  the  forecast. 

4.  Multiplying,  adding,  or  determining  trends  among  predietors  added  minimal 
accuracy. 
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TABLE  2.  List  of  Predictands.  Vislocky  and  Fritsch  (1997)  produced  forecasts 
using  the  surfaee  observational  based  system  for  these  signifieant  aviation 
operational  thresholds.  Adapted  from  Vislocky  and  Fritsch  (1997). 


Variable 

Binary  Thresholds 

Ceiling  Height 

<  500  ft 

<  1000  ft 

<  3000  ft 

<  6500  ft 

Visibility 

<  1  mi 

<  3  mi 

<  5  mi 

<7  mi 

5.  The  optimal  number  of  stations  to  use  for  1,  3,  and  6-hr  forecasts  are  10,  25, 
and  40  respectively. 

Visloeky  and  Fritseh  (1997)  determined  that  33  surfaee  parameters  were 
operationally  signifieant  in  forecasting  fog  dissipation  (Table  3).  They  also  eoneluded 
that  the  addition  of  climatologieal  faetors  relating  to  incoming  solar  radiation  produced  a 
more  accurate  forecast. 


TABLE  3.  Surface  Meteorological  Parameters  Used  as  Predictors.  These  33  conditions 
were  determined  to  be  the  most  significant  predictors  of  fog  dissipation.  Adapted  from 
Vislocky  and  Fritsch  (1997). 


Variable 

Binary  Threshold 

Opaque  Cloud  Amount 

>  1/10 

>5/10 

>9/10 

Total  Cloud  Cover 
Precipitation  Occurrence 

Clear 

Yes 

Scattered 

Broken 

Overcast 

Obscured 

Ceiling  Height 

<  200  ft 

<  500  ft 

<  1000  ft 

<  3000  ft 

<  6500  ft 

<  12000  ft 

Visibility 

<  0.5  mi 

<  1.0  mi 

<3.0  mi 

<  5.0  mi 

<  7.0  mi 

<  10  mi 

Wind  Direction 

Wind  Direction  (cont.) 
Dew  Point 

Dewpoint  Depression 

Sea  Level  Pressure 

0° 

<270° 

<45° 

<315° 

<90° 

<  135° 

<  180° 

<235° 

9 


Vislocky  and  Fritsch  (1997)  used  a  least  squares  linear  regression  model  to 
develop  their  predietive  equations.  In  this  method,  eaeh  predietand  is  the  result  of  a 
eombination  of  predictors.  A  forward  stepwise  screening  algorithm  was  used  to  select 
the  best  available  predictors  (Vislocky  and  Fritsch  1997).  One  predictor  is  added  at  a 
time  in  order  to  gauge  the  effect  each  has  on  the  predietand.  This  method  allows  useless 
or  redundant  predictors  to  be  discarded  (Vislocky  and  Fritsch  1997). 

The  results  of  Vislocky  and  Fritsch’s  (1997)  findings  can  be  summarized; 

1.  Observations-based  methods  on  average  produced  12%  (with  a  5-20%  range) 
better  results  than  conditional  climatology  and  4%  better  than  MOS. 

2.  Any  short-range  forecast  should  rely  more  heavily  on  observations  than 
numerical  models. 

3.  At  the  6-hour  point,  the  observations  and  MOS  predictands  were  similar, 
which  verified  that  the  6-hr  point  is  the  time  for  numerical  models  to  take 
over. 

Two  major  advantages  of  an  observations-based  system  are  that  it  is  not  tied  to 
any  numerical  model  and  it  can  be  run  in  a  matter  of  seconds  on  any  computer  as  soon  as 
new  data  becomes  available  (Vislocky  and  Fritsch  1997). 

Suggested  improvements  to  the  observations-based  statistical  method  include 
incorporation  of  radar  data,  satellite  imagery,  and  upper  air  observations  into  the 
equations  (Vislocky  and  Fritsch  1997).  Hilliker  and  Fritsch  (1999)  added  upper  air 
observations  to  the  surface-based  observations  system  in  a  study  at  the  San  Francisco 
International  Airport.  They  used  upper  air  data  from  Oakland  since  it  is  within  the  area 
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of  interest  and  developed  additional  upper  air  predictors  (Table  4).  Static  stability  was 
added  as  a  predictor  in  addition  to  basic  upper  air  parameters. 

Hilliker  and  Fritsch  (1999)  used  a  predictand  of  3000  ft  ceiling  height  and  surface 
predictors  similar  to  Vislocky  and  Fritsch  (1997)  with  modifications  for  local  effects. 
Hilliker  and  Fritsch  (1999)  employed  a  logarithmic  regression  model  with  the  addition  of 
upper  air  parameters.  Hilliker  and  Fritsch’s  (1999)  research  showed  that: 

1.  Including  upper  air  predictors  produced  0-3%  improvements  over  the  strictly 
surface  based  systems  for  0-3  hr  forecasts. 

2.  The  longer  the  forecast  time,  the  less  impact  the  upper  air  observations  had  on 
the  forecast. 

3.  Improvements  of  up  to  32%  were  shown  over  MOS  forecasts. 

Hilliker  and  Fritsch  (1999)  proved  that  even  in  a  region  with  limited  observations, 
observations-based  systems  are  superior  to  conditional  climatology,  persistence,  and 
MOS  forecasts  for  short-term  forecasts.  Furthermore,  including  upper  air  parameters, 
such  as  static  stability,  increased  the  accuracy  of  the  forecasts.  Two  additional 
parameters  that  could  improve  the  forecast  are  the  inversion  height  and  the  thickness  of 
the  cloud  layer  (Hilliker  and  Fritsch  1999). 

Vislocky  and  Fritsch  (1997)  and  Hilliker  and  Fritsch  (1999)  showed  success  in 
producing  dissipation  forecasts  using  observations-based  networks.  The  next  section 
examines  the  effectiveness  of  evaluating  satellite  imagery  to  forecast  dissipation. 
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TABLE  4.  Upper  Air  Meteorological  Parameters  Used  as  Predictors.  Parameters  were 
considered  for  8  pressure  levels  between  1000  mb  and  500  mb,  with  the  exception  of 
static  stability,  which  was  evaluated  in  three  layers  between  1000  mb-850  mb.  Adapted 
from  Hilliker  and  Fritsch  (1999). 


Parameter 

Binary  Threshold 

Height 

Temperature 

Relative  Humidity 

Wind  Direetion 

Wind  Direetion  (cont.) 
Wind  Speed 

Static  Stability  (d0/dz) 

>  30% 

23°  to  <68° 

203°  to  <  248° 

>  50%  >  70% 

68°  to  <113°  113°to<158° 

248°  to  <293°  293°  to  <338° 

>  90% 

158°  to  <203° 

338°  to  <23° 

2.4  Satellite  Imagery 

Another  method  in  predicting  the  dissipation  of  fog  is  the  use  of  satellite  imagery. 
Gurka  (1974,  1978)  showed  that  high  resolution  visible  imagery  can  be  used  to  predict 
fog  dissipation  because  of  the  strong  correlation  between  fog  brightness  and  duration. 
Dissipation  occurs  first  on  the  outer  edges  of  the  fog  (Gurka  1974),  where  the  ground  is 
heated  more  intensely  along  the  boundaries  than  in  the  interior  (Gustafson  and 
Wasserman  1976).  Anthis  and  Cracknell  (1999)  concluded  that  the  dissipation  of  fog 
proceeds  inward  with  time  and  that  it  is  dissipated  according  to  layer  thickness. 

Gurka’ s  (1978)  work  provided  significant  findings,  concluding  that  enhanced 
visible  satellite  imagery  could  produce  effective  short-range  forecasts  of  fog  dissipation, 
especially  for  aviation.  Gurka  (1978)  also  realized  that  computers  would  be  able  to 
ingest  and  process  this  data  quickly  enough  to  produce  real  time  short-range  predictions. 
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Since  all  previously  mentioned  methods  had  suceess  in  predicting  dissipation 
time,  the  next  seetion  examines  the  possibility  of  eombining  parameters  into  a  multi- 
souree  method. 

2.5  Multi-Source  Methods 

Reudenbach  and  Bendix  (1998)  made  progress  using  multi-souree  data  to  prediet 
dissipation  in  Germany.  Their  work  used  four  separate  types  of  data  for  model 
development:  thermodynamic  equations,  NOAA  Advaneed  Very  High  Resolution 
Radiometer  (AVHRR)  imagery,  terrain  models,  and  surfaee  and  upper  air  meteorologieal 
observations.  Reudenbaeh  and  Bendix  (1998)  foeused  on  multi-souree  data  beeause 
single  souree  methods  are  not  able  to  aeeount  for  all  aspeets  of  dissipation.  Of  note  is 
that  although  Gurka  (1974)  presented  a  method  for  noweasting  fog  clearanee,  the  method 
was  never  standardized,  nor  proven  100%  reliable.  Reudenbaeh  and  Bendix  (1998)  also 
noted  that  numerieal  models  eannot  always  aeeount  for  loeal  effeets  sueh  as  dissipation 
due  to  advective  proeesses.  They  suceessfully  developed  a  straightforward  model  for 
foreeasting  elearanee  in  1km  by  1km  areas. 

This  ehapter  summarizes  multiple  methods  employed  to  aeeurately  prediet  the 
dissipation  of  fog.  Visloeky  and  Fritseh  (1997)  and  Hilliker  and  Fritseh  (1999)  showed 
that  an  observations-based  method  produees  superior  forecasts,  while  standards  sueh  as 
MOS  guidanee,  persistenee  and  eonditional  elimatology  are  still  regularly  employed. 
Gurka  (1974,  1978)  and  Anthis  and  Craeknell  (1998,  1999)  demonstrated  that  visible 
satellite  imagery  eould  be  used  to  prediet  duration  of  fog  events.  Reudenbaeh  and 
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Bendix  (1998)  demonstrated  the  effeetiveness  of  a  multi-souree  data  system  on  the 
aeeurate  predietion  of  fog  dissipation  times 

Eaeh  method  displayed  some  degree  of  sueeess  in  predieting  fog  dissipation; 
however  this  researeh  foeuses  on  adapting  a  statistieally-based  method  to  foreeast 
dissipation.  There  are  several  reasons  for  this.  First,  as  shown  by  Visloeky  and  Fritseh 
(1997)  this  methodology  is  superior  to  model  output  and  eonditional  elimatology. 
Seeond,  satellite  imagery  by  itself  has  not  proven  to  be  a  single  effeetive  tool  for 
predieting  dissipation.  A  multi-souree  method,  that  ineludes  an  observations-based 
network  and  satellite  imagery,  eould  produee  a  superior  foreeast  and  should  be  evaluated 
in  the  future.  The  next  ehapter  details  the  proeess  of  developing  and  implementing  the 
statistieally-based  system. 
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III.  Methodology 


3.1  Overview 

The  literature  review  deseribes  several  methods  for  predicting  fog  dissipation  that 
prove  effective  in  all  geographic  areas.  Vislocky  and  Fritsch  (1997)  showed 
observations-based  systems  can  be  successfully  employed  along  the  east  coast  of  the  US, 
while  Hilliker  and  Fritsch  (1999)  showed  that  the  system  is  effective  for  the  west  coast. 
Since  observations-based  networks  are  superior  to  conditional  climatology  and  model 
output  (Vislocky  and  Fritsch  1997),  this  method  is  applied  to  the  airfields  in  the  15  OWS 
to  produce  a  superior  fog  dissipation  forecast  tool. 

This  research  follows  Vislocky  and  Fritsch  (1997)  and  Hilliker  and  Fritsch’s 
(1999)  work  into  the  use  of  a  statistically-based  observation  network  to  produce  accurate 
short-term  dissipation  forecasts.  Although  this  method  is  a  proven  success  along  the  east 
coast,  changes  are  made  to  fit  the  needs  of  the  AF.  Specifically,  a  different  set  of 
predictands  is  required  in  order  to  meet  the  need  of  the  15  OWS,  as  Vislocky  and  Fritch’s 
(1997)  original  work  failed  to  address  the  critical  threshold  of  200  ft  ceilings  and  0.5  mi 
visibility. 

The  scope  of  the  problem  for  the  entire  15  OWS  AOR  may  seem  large,  however, 
it  is  not  insurmountable  and  is  solved  piecewise.  The  most  critical  problem,  200  ft 
ceilings  and  0.5  mile  visibility  on  the  east  coast,  is  addressed  in  this  research.  Once  the 
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technique  and  processes  are  established,  equations  tested,  and  methodology  set,  other 
airfields  and  predictands  can  be  examined  in  the  future. 

The  remaining  portion  of  this  chapter  is  dedicated  to  an  explanation  on  the 
procedures  and  processes  employed  to  develop  the  predictive  equations.  Section  3.2 
examines  the  surface,  upper  air,  and  climatological  data  sets  used  to  develop  the 
equations,  while  section  3.3  details  the  statistical  approach  taken  in  this  research. 

3.2  Data 

3.2.1  Surface  and  Upper  Air  Data.  This  research  is  based  on  a  five-year  data  set  (1998- 
2002),  consisting  of  standard  surface  observations  in  METAR  format  and  raw  upper  air 
data  in  formatted  ASCII  text,  provided  by  the  Air  Force  Combat  Climatology  Center 
(AFCCC)  in  Asheville,  NC.  The  surface  data  consist  of  both  standard  hourly  and  special 
observations  for  stations  along  the  east  coast  of  the  US.  The  upper  air  data  consist  of 
both  mandatory  and  significant  levels  and  are  available  twice  a  day,  OOZ  and  12Z. 

The  large  data  set  is  required  for  two  reasons.  First,  due  to  the  relatively  rare 
occurrence  of  the  extreme  conditions  (200  ft  ceilings/  0.5  mile  visibility),  a  large  data  set 
ensures  a  statistically  significant  number  of  occurrences.  Second,  the  data  are  broken 
into  two  subsets  in  order  to  develop  the  predictive  equations  with  one  set  and 
independently  verify  the  results  with  another. 

The  implementation  of  the  Automated  Surface  Observing  System  (ASOS) 
program  by  the  NWS  in  the  middle  to  late  1990’s  had  a  significant  effect  on  this  research. 
The  advantages  of  the  automated  systems  are  significant.  The  density  of  the  surface 
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observation  network  in  the  US  inereased  signifieantly  both  spatially  and  temporally. 
Increased  spatial  coverage  results  in  smaller  surface-based  networks,  which  in  turn  are 
more  representative  of  the  conditions  at  the  airfield  of  interest.  Increased  temporal 
resolution  occurs  as  many  part  time  observing  stations  (which  were  excluded  from 
Vislocky  and  Fritsch’s  (1997)  initial  study)  were  transitioned  into  full  time  observing 
stations  in  the  late  90’s.  This  positively  influences  this  study  in  two  ways.  First,  stations 
which  were  previously  ignored,  could  now  be  evaluated,  once  again  increasing  the  spatial 
density.  Second,  more  data  are  available  for  all  hours  of  the  day,  resulting  in  evaluation 
of  more  cases  of  fog  occurrence  and  dissipation.  While  these  advantages  assist  in  this 
work,  problems  also  arise. 

First,  the  system  was  not  implemented  instantaneously  at  each  location.  That  is, 
there  is  a  significant  difference  between  the  times  in  which  each  airfield’s  system  came 
online.  For  this  reason,  stations  fitted  with  ASOS’s  after  1998  are  omitted  from  this 
study,  even  if  they  fell  within  the  observations-based  network,  due  to  their  lower  data 
density.  Second,  the  initial  period  of  data  received  from  the  ASOS  is  often  missing 
weather  variables.  These  are  two  of  the  main  factors  for  the  development  and 
implementation  of  missing  data  schemes. 

3.2.2  Missing  Surface  Data.  Missing  data  are  broken  up  into  three  different  categories, 
each  of  which  is  handled  in  a  separate  way.  The  goal  of  the  missing  data  algorithms  is 
two-fold.  First,  a  scheme  is  developed  to  maximize  the  number  of  data  points  available 
for  the  study.  Second,  the  replacement  data  needs  to  reflect  the  current  meteorological 
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conditions  at  each  station  as  accurately  as  possible.  For  this  reason  a  three  step  missing 
data  approaeh  is  implemented. 

Step  1 :  Missing  eeiling  data  due  to  the  oecurrenee  of  fog.  One  of  the  most 
eommon  oeeurrenees  of  missing  data  is  the  laek  of  sky  eondition  and/or  eeiling  height 
when  fog  is  reported.  For  this  reason  eeiling,  values  are  represented  aeeording  to  Table 
5,  data  are  missing  with  redueed  visibility  due  to  the  oeeurrence  of  fog. 

Step  2:  Missing  entire  observation.  If  an  entire  surfaee  observation  is  missing,  a 
nearest  reliable  neighbor  approach  is  used.  A  reliable  neighbor  is  defined  as  a  full  time 
observing  station  throughout  the  entire  data  set.  The  nearest  neighbor  approaeh  is  used 
beeause  it  is  the  most  physically  significant  replacement  proeess  available.  Fog  generally 
oeeurs  in  the  eold  season  and  is  eonsidered  more  of  a  synoptie  seale  rather  than 
mesoseale  event.  For  this  reason  surfaee  weather  eonditions  at  nearby  loeations  are 
typieally  similar  enough  (espeeially  along  the  data  dense  east  eoast)  to  make  this  a  sound 
meteorologieal  argument.  The  exeeption  is  redueed  visibility  or  eeiling  conditions  due  to 
loealized  eonveetive  events  or  loealized  terrain  infiuenees,  whieh  are  assumed  to  oeeur 
infrequently  enough  to  not  have  a  signifieant  impact  on  this  study. 


TABLES.  Binary  Ceiling  Replacement  Values.  The 
replaeement  values  used  for  ceilings  when  the  sky  condition  is 
missing  from  an  observation  and  fog  is  reported. 


Observed  Visibility 

Binary  Ceiling  Value 

Vis  <  0.5  mi 

Cig  <  200  ft 

Vis  <  1.5  mi 

Cig  <  1000  ft 

Vis  <  3  mi 

Cig  <  3000  ft 
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Step  3:  Missing  individual  weather  data  from  an  observation.  This  type  of 
missing  data  oecurs  when  single  or  multiple  elements  from  a  surface  observation  are 
missing.  Data  are  replaced  once  again  using  the  nearest  neighbor  approach,  but  this  time 
only  replacing  the  individual  element. 

The  missing  data  schemes  are  only  applied  to  predictors  and  never  to  the 
predictand.  If  the  predictand  is  missing  at  any  time,  that  occurrence  is  not  included  in  the 
study  for  either  model  building  or  verification  purposes.  After  the  surface  data  is 
replaced  missing  upper  air  data  algorithms  are  implemented. 

3.2.3  Missing  Upper  Air  Data.  Missing  upper  air  data  are  handled  quite  differently  than 
missing  surface  data  due  to  the  relative  sparseness  of  observations  both  spatially  and 
temporally.  Missing  data  are  replaced  using  data  from  the  upper  air  observation  24  hours 
earlier.  Therefore,  OOZ  upper  air  observations  are  replaced  with  OOZ  data  the  day  prior. 

If  the  data  are  unavailable  24  hours  ago,  the  event  is  not  included  in  the  study.  Data  from 
the  previous  24  hours  are  used  in  order  to  best  represent  atmospheric  conditions.  Since 
observations  are  only  taken  twice  a  day,  the  OOZ  data  are  not  typically  representative  of 
the  12Z  observations  and  visa  versa. 

This  technique  is  applied  as  an  accurate  representation  based  on  two  factors. 

First,  the  synoptic  situation  typically  associated  with  fog  events  is  stable  and  not  rapidly 
changing.  Second,  replacement  algorithms  are  used  infrequently  enough  in  most  cases  to 
not  have  a  significant  effect  on  this  research.  The  Dover  AFB  and  McGuire  AFB  data 
sets  each  have  3.7%  missing  data  rates.  The  Andrews  AFB  data  set  has  a  12%  missing 
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data  rate,  this  case  is  a  concern  heading  into  equation  development.  Next  examined  is 
the  processing  of  the  data. 

3.2.4  Data  Assimilation  and  Manipulation.  The  first  step  in  this  procedure  is  setting  up 
the  surface  networks  for  each  of  the  airfields  for  each  of  the  six  hour  increments. 

Vislocky  and  Fritsch  (1997)  concluded  that  the  10,  25,  and  40  closest  observing  stations 
are  the  optimal  number  of  stations  to  be  used  in  1,  3,  and  6  hour  forecasts.  The  optimal 
number  of  stations  is  interpolated  as  19,  30,  and  35  for  the  2,  4,  and  5  hour  forecasts. 
Observational  networks  are  built  independently  for  each  forecast  location  and  time.  The 
surface  observing  stations  for  the  observational  networks  included  in  this  study  are  for 
Dover  AFB  (KDOV),  McGuire  AFB  (KWRI),  and  Andrews  AFB  (KADW)  and  are 
detailed  in  the  Appendix. 

Data  processing  is  the  main  challenge  of  this  exercise  and  is  summarized  in 
Figure  2.  The  first  step  is  an  evaluation  of  the  airfield  of  interest  for  occurrence  of  the 
predictand  value  (e.g.  ceiling  less  than  200  ft).  If  an  observation  has  an  occurrence  below 
one  of  the  predictand  thresholds,  it  is  extracted.  The  same  observation  may  meet  multiple 
criteria  (i.e.  ceiling  below  200  ft  and  visibility  less  than  0.5  mi).  In  this  case,  the 
observation  is  used  for  all  applicable  predictand  criteria.  This  experiment  is  set  up  to 
evaluate  every  case  of  fog  that  meets  the  criteria  regardless  of  time  or  season. 

The  next  step  is  an  evaluation  of  network  observations.  The  date  time  group  of 
each  of  the  40  stations  in  the  observations  network  is  compared  to  the  main  station’s  date 
time  group.  Observations  within  an  hour  are  extracted  to  a  separate  file.  If  no 
observation  is  available  within  an  hour,  it  is  left  for  missing  data  algorithms  to  replace. 
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The  same  process  is  then  applied  to  the  upper  air  and  climatological  data.  With  the 
extraction  of  the  relevant  data,  the  statistical  approach  is  now  reviewed. 


3.3  Statistical  Approach 


3.3.1  Predictands.  After  the  network  is  constructed  new  predictands  are  evaluated.  As 
noted  earlier,  the  predictands  offered  by  Vislocky  and  Fritsch  (1997)  do  not  fit  the  needs 


FIG  2.  Flowchart  for  Data  Manipulation. 
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of  the  15  OWS,  where  the  eritical  thresholds  (200  ft  eeilings  and  0.5  mi  visibility)  are  not 
addressed.  This  work  evaluates  predietands  of  200  ft  eeiling  and  0.5  mile  visibility  at 
three  airfields  along  the  eastern  seaboard  of  the  US:  MeGuire  AFB,  Dover  AFB,  and 
Andrews  AFB.  Predietands  are  assigned  a  binary  value  of  ‘  1  ’  for  oeeurrenee  of  an  event 
and  ‘0’  for  non-oeeurrenee  of  event.  Predietive  equations  are  developed  for  each 
individual  predictand  at  each  airfield  at  hourly  increments  spanning  one  to  six  hours 
beyond  the  observation  of  fog  occurrence.  The  predictors  that  influence  the  predictand 
forecast  are  now  examined. 

3.3.2  Predictors.  This  investigation  deviates  from  Vislocky  and  Fritsch’s  (1997)  work  in 
terms  of  surface  predictors.  The  main  discrepancy  is  opaque  cloud  cover,  which  is 
omitted  from  this  study.  The  reason  for  this  is  twofold.  First,  Vislocky  and  Fritsch’s 
(1997)  work  used  an  earlier  data  set  consisting  of  surface  observations  in  SA  (surface 
airways)  format.  These  observations  include  more  ceiling  information  than  is  currently 
reported  in  METAR  observations  (namely  the  presence  or  absence  of  a  ceiling). 
Furthermore,  SA  ceilings  were  more  subjective,  that  is,  even  if  8/10  of  the  sky  were 
covered  by  clouds,  an  observer  could  determine  it  to  not  be  a  ceiling  due  to  the  lack  of 
opacity.  The  second  reason  is  inconsistency  in  remark  data,  especially  reporting  the  eight 
and  nine  cloud  groups  which  detail  cloud  amount  and  type.  Since  these  are  not  present 
on  a  majority  of  observations,  it  is  impossible  to  determine  with  the  current  available  data 
the  opacity  of  the  sky. 

A  second  difference  in  this  research  from  Vislocky  and  Fritsch’s  (1997)  is  the  use 
of  altimeter  setting  rather  than  sea  level  pressure.  The  reason  for  this  is  that  altimeter  is 
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reported  in  all  observations  and  sea  level  pressure  is  not  always  reported  in  the  remarks 
of  the  data  set.  Sinee  all  stations  are  within  a  eouple  of  hundred  meters  of  sea  level,  the 
approximation  of  altimeter  setting  for  sea  level  pressure  is  valid.  Binary  predietors  are 
eoded  with  a  “dummy”  value  of  ‘0’  representing  non-oeeurrenee  or  ‘  1  ’  representing 
oeeurrenee,  and  non-binary  predietors  maintain  their  observed  values.  Eaeh  set  of  30 
predietors  (Table  6)  is  evaluated  at  eaeh  surfaee  observing  station  within  the  observing 
network  (the  number  of  stations  in  eaeh  observations  network  varies  with  time). 


TABLE  6.  Surfaee  Predietors  for  the  15  OWS  AOR.  The  surfaee  meteorologieal 
predietors  used  to  ereate  the  observations-based  foreeasts  for  the  15  OWS  AOR.  All 
predietors  are  evaluated  at  eaeh  airfield  in  the  speeifie  network  (whieh  varies  depending 
on  the  time  of  the  foreeast). 


Variable 

Designator 

Occurrence 

Binary 

Non  occurrence 

Non-Binary 

Sky  Condition:  Clear 

SKY  (SKC) 

I 

0 

Sky  Condition:  Few/Scattered 

SKY  (FEW/SCT) 

I 

0 

Sky  Condition:  Broken 

SKY  (BKN) 

I 

0 

Sky  Condition:  Overcast 

SKY  (OVC) 

I 

0 

Sky  Condition:  Obscured 

SKY  (OBSCURED) 

I 

0 

Precipitation 

PRECIP 

I 

0 

Ceiling  Height  <  200  ft 

CIG  HT  <200  FT 

I 

0 

Ceiling  Height  <  500  ft 

CIG  HT  <500  FT 

I 

0 

Ceiling  Height  <  1000  ft 

CIG  HT  <1000  FT 

I 

0 

Ceiling  Height  <  3000  ft 

CIG  HT  <3000  FT 

I 

0 

Ceiling  Height  <  6500  ft 

CIG  HT  <6500  FT 

I 

0 

Ceiling  Height  <  12000  ft 

CIG  HT  <12000  FT 

I 

0 

Visibility  <  0.5  mi 

VIS  <  800  M 

I 

0 

Visibility  <  I  mi 

VIS  <  1600  M 

I 

0 

Visibility  <  3  mi 

VIS  <  4800  M 

I 

0 

Visibility  <  5  mi 

VIS  <  8000  M 

I 

0 

Visibility  <  7  mi 

VIS  <  9999  M 

I 

0 

Visibility  <  10  mi 

VIS  >=  9999  M 

I 

0 

Wind  Direction:  Variable 

WIND  DIR  (VRB) 

I 

0 

Wind  Direction  <  45 

WIND  DIR  (<45) 

I 

0 

Wind  Direction  <  90 

WIND  DIR  (<90) 

I 

0 

Wind  Direction  <135 

WIND  DIR  (<135) 

I 

0 

Wind  Direction  <  180 

WIND  DIR  (<I80) 

I 

0 

Wind  Direction  <  225 

WIND  DIR  (<225) 

I 

0 

Wind  Direction  <  270 

WIND  DIR  (<270) 

I 

0 

Wind  Direction  <315 

WIND  DIR  (<3I5) 

I 

0 

Wind  Direction  <  360 

WIND  DIR  (<360) 

I 

0 

Dewpoint 

DEWPOINT 

X 

Dewpoint  Depression 

DEWPOINT  DEPRESSION 

X 

Altimeter  Setting 

ALT 

X 
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The  upper  air  predietors  follow  the  work  of  Hilliker  and  Fritseh  (1999)  and  are 
available  in  Table  7.  Eaeh  predietor  is  evaluated  at  six  levels,  surfaee,  1000  mb,  925  mb, 
850  mb,  700  mb,  and  500  mb,  with  the  exeeption  of  statie  stability,  whieh  is  evaluated  in 
three  layers:  1000  mb-850  mb,  1000  mb-925  mb,  and  925  mb-850  mb.  The  binary 
predictors  are  coded  in  the  same  manner  as  the  surface  predictors.  Upper  air  data  from 
OOZ  are  used  for  fog  occurrence  times  of  03Z-15Z  (12Z  data  are  used  for  the  hours  15Z- 
03Z)  in  order  to  represent  lag  times  in  data  availability,  ensuring  this  method  can  be 
implemented  operationally.  The  upper  air  sounding  that  is  closest  to  each  airfield  of 
interest  is  used  (Table  8)  and  is  considered  to  be  representative  of  the  upper  air  conditions 
for  the  entire  surface  observation  network. 

Numerous  climatological  predictors  are  available  for  selection  by  the 
predictive  equations.  First  among  these  are  climatological  frequencies  of  occurrence  of 


TABLE  7.  Upper  Air  Predictors  for  the  15  OWS  AOR.  The  upper  air  meteorological 
predictors  used  to  create  the  observations-based  forecasts  for  the  15  OWS  AOR.  All 
predictors  are  evaluated  at  each  pressure  with  the  exception  of  static  stability,  which  is 
evaluated  in  the  three  layers. 


Variable 

Designator 

Occurrence 

Binary 

Non  occurrence 

Non-Binary 

Height 

HT 

X 

Temperature 

T 

X 

Relative  Humidity  >  30% 

RH  >  30% 

1 

0 

Relative  Humidity  >  50% 

RH  >  50% 

1 

0 

Relative  Humidity  >  70% 

RH  >  70% 

1 

0 

Relative  Humidity  >  90% 

RH  >  90% 

1 

0 

Wind  Direction:  23°  to  <  68° 

WINDN 

1 

0 

Wind  Direction:  68°  to  <  1 13° 

WIND  NE 

1 

0 

Wind  Direction:  1 13°  to  <  158° 

WINDE 

1 

0 

Wind  Direction:  158°  to  <  203° 

WIND  SE 

1 

0 

Wind  Direction:  203°  to  <  248° 

WINDS 

1 

0 

Wind  Direction:  248°  to  <  293° 

WIND  SW 

1 

0 

Wind  Direction:  293°  to  <  338° 

WIND  W 

1 

0 

Wind  Direction:  338°  to  <  23° 

WIND  NW 

1 

0 

Wind  Speed 

WIND  SPEED 

X 

Static  Stabililty  (d9/dz):  (1000-925  mb) 

STATIC  STABILITY  (1000-925  MB) 

X 

Static  Stabililty  (dO/dz):  (1000-850  mb) 

STATIC  STABILITY  (1000-850  MB) 

X 

Static  Stabililty  (dO/dz):  (925-850  mb) 

STATIC  STABILITY  (925-850  MB) 

X 
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TABLE  8.  Upper  Air  Stations  in  the  15  OWS  AOR.  The  three  upper 
air  stations  used  in  equation  development  along  with  the  airfield  they 
represent. 


Airfield 

ICAO 

Upper  Air  Station 

ICAO 

McGuire  AFB  (NJ) 

KWRI 

Brookhaven  (NY) 

KOKX 

Andrews  AFB  (VA) 

KADW 

Washington  DC 

KIAD 

Dover  AFB  (DE) 

KDOV 

Wallops  Island  (MD) 

KWAL 

eeiling  and  visibility  threshold  at  the  airfield  of  interest.  Predictor  values  including  200 
ft/0.25  mi,  500  ft/0.5  mi,  1000  ft/1.5  mi  are  extracted  from  the  Operational  Climatic  Data 
Summary  tables  available  through  AFCCC.  Sine  and  cosine  of  the  day  of  the  year  are 
also  included  as  predictors  of  solar  radiation  effects  per  Vislocky  and  Fritsch  (1997). 

A  better  understanding  of  the  effect  of  solar  radiation  on  the  dissipation  process  is 
also  included.  Following  the  work  of  Campbell  (1977),  the  maximum  amount  of 
incoming  solar  radiation  reaching  the  surface  of  the  earth  for  the  time  of  year  is 
calculated.  These  calculations  do  not  take  into  account  clouds  above  the  fog  and  stratus; 
however,  for  this  study  it  is  assumed  that  all  of  this  radiation  is  reaching  the  top  of  the 
cloud/fog  layer.  All  of  the  climatological  predictors  are  non-binary,  and  a  summary  is 
available  in  Table  9.  Now  that  predictors  and  predictands  are  explained  the  statistical 
methods  behind  equation  development  is  examined. 

3.3.3  Equation  Development.  A  least  squares  multiple  linear  regression  model  is  used  to 
develop  six  predictive  equations  for  each  of  the  two  predictands  at  each  airfield.  The 
first  assumption  made  in  order  to  use  linear  regression  is  that  errors  are  uncorrelated 
random  variables  with  constant  variance  and  zero  for  a  mean  (defined  as 
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TABLE  9.  Climatological  Predictors  for  the  15  OWS  AOR.  The  climatological 
predictors  used  to  create  the  observations-based  forecasts  for  the  15  OWS  AOR.  All 
predictors  are  evaluated  at  only  the  airfield  of  interest  (either  Dover  AFB,  McGuire  AFB, 
or  Andrews  AFB). 


Variable 

Designator 

Occurrence 

Binary 

Non  occurrence 

Non-Binary 

Frequency  of  Occurrence:  200  ft/0.25  mi 

CLIMO  <  200/0#25 

X 

Frequency  of  Occurrence:  500  ft/0.5  mi 

CLIMO  <  500/0#5 

X 

Frequency  of  Occurrence:  1000  ft/1.5  mi 

CLIMO  <  1000/1#5 

X 

Sine  of  the  Day 

SIN  DAY 

X 

Cosine  of  the  Day 

COS  DAY 

X 

Solar  Radiation 

SOLAR  RADIATION 

X 

homoscedastieity)  (Montgomery  and  Runger  2003).  Analysis  of  the  residuals  from  eaeh 
linear  model  can  prove  these  assumptions  valid.  The  residual  of  a  regression  model  is 
defined; 

ei=yi-yi  1=1,2...,  n  (1) 

Where  yt  is  the  observed  value  and  y,  is  the  fitted  value  obtained  from  the  linear 
regression  model.  A  randomly  distributed  residual  seatter  plot,  e.g.  showing  no  pattern  or 
trend,  verifies  the  above  assumptions. 

Figure  3  shows  sample  residual  scatter  plots  for  eaeh  of  the  three  data  sets.  The 
above  eharacteristie,  a  random  seattering  of  points  around  the  zero  line,  is  present  in  the 
diagrams  below.  This  leads  to  the  conclusion  that  the  models  have  a  mean  of  zero  and 
eonstant  variance.  The  above  test  is  performed  on  eaeh  of  the  36  models  individually. 

The  seeond  assumption  required  to  use  linear  regression  is  that  the  errors  are 
normally  distributed  (Montgomery  &  Runger  2003).  First,  the  Central  Fimit  Theorem 
states  that  as  a  sample  size  gets  larger  (typically  assumed  to  be  greater  than  30),  the 
distribution  tends  towards  normal  (Montgomery  &  Runger  2003).  A  more  precise  way  to 
eheck  the  distribution  is  by  examining  a  normal  probability  plot  of  the  residuals  for  each 
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FIG.  3.  Residual  Plots  from  Data  Sets.  Residual  plots  from  the  KDOV  Cig  LT  200  ft 
(T+2)  (a),  KADW  Cig  LT  200  ft  (T+2)  (b),  and  KWRI  Vis  LT  0.5  mi  (T+2)  (c)  data  sets. 
These  typical  patterns  show  residual  plots  scattered  and  randomly  distributed  around  the 
zero  line,  confirming  the  assumptions  necessary  to  use  linear  regression. 


of  the  36  models.  To  assume  normality,  the  residuals  should  fall  along  a  straight  line. 
Figure  4  shows  two  examples  of  normal  probability  plots.  The  first  is  from  the  Andrews 
AFB  ceiling  less  than  200  ft  forecast  at  one-hour  (Fig.  4a).  Notice  the  plot  is  in  almost  a 
normal  distribution  of  residuals.  The  second  example  is  from  the  Dover  AFB  ceiling  less 
than  200  ft  forecast  at  two-hours  (Fig.  4b).  Notice  the  plot  falls  essentially  along  a 
straight  line,  allowing  the  conclusion  of  normality  (Montgomery  &  Runger  2003). 

Table  10  lists  the  36  derived  models  for  this  work  and  details  the  statistical  tests  run 
for  each.  Notice  that  10  of  the  36  models  have  errors  which  are  not  normally  distributed 
(bold,  column  2,  Table  10).  Six  of  the  non-normal  distributions  are  found  in  the  Andrews 
AFB  data  set.  Despite  these  discouraging  indicators,  regression  models  are  nonetheless 
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a) 

Normal  Quantile  Plot 


b) 


Normal  Quantile  Plot 


FIG.  4.  Normal  Probability  Plots.  Plots  for  KADW  Cig  LT  200  ft  (T+1)  (a)  and  KDOV 
Cig  LT  200  ft  (T+2)  (b).  Note  the  non-linear  shape  of  the  first  model,  leading  to  the 
conclusion  of  a  lack  of  normalcy.  Panel  (b)  exhibits  a  linear  pattern,  indicating  a  normal 
distribution  of  the  residuals,  which  is  a  necessary  assumption  for  linear  regression, 
stair  step  pattern,  not  close  to  resembling  a  straight  line.  This  is  an  example  of  a  non- 


developed  and  analysis  of  effects  of  the  non-normal  distributions  is  found  in  Chapter  IV. 
A  lack  of  homoscedasticity  (e.g.  heteroscedasticity)  is  present  in  only  three  of  the  models 
(bold,  column  3,  Table  10).  However,  as  Neter  et  al.  (1990)  point  out,  if 
heteroscedasticity  is  inherent,  the  results  obtained  are  still  unbiased,  consistent 
estimators,  but  no  longer  possess  minimum  variance.  Since  the  KWRI  Vis  LT  0.5  mi 
(T+1)  and  KADW  Cig  LT  200  ft  (T+1)  possess  neither  homoscedacity  nor  normal 
distributions,  these  data  are  of  particular  concern  going  into  model  development. 

Multiple  linear  regression  examines  the  relationship  between  two  or  more  variables 
(Montgomery  and  Runger  2003)  with  each  equation  taking  the  form; 

Y=|3o+(3iXi+(32X2+...+(3nXn+E  (2) 
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TABLE  10.  Statistical  Tests.  Definitions  of  the  36  models  derived  in 
this  study  along  with  the  results  of  the  statistieal  tests  for  normality  and 
heteroseedastieity  needed  to  eonfirm  assumptions  neeessary  to  use 
linear  regression. 


Forecast 

Normal 

Homoscedasticity 

KDOV  Cig  LT  200  ft  (T+1) 

Yes 

Yes 

KDOV  Cig  LT  200  ft  (T+2) 

Yes 

Yes 

KDOV  Cig  LT  200  ft  (T+3) 

Yes 

Yes 

KDOV  Cig  LT  200  ft  (T+4) 

Yes 

Yes 

KDOV  Cig  LT  200  ft  (T+5) 

Yes 

Yes 

KDOV  Cig  LT  200  ft  (T+6) 

Yes 

Yes 

KDOV  Vis  LT  0.5  mi  (T+1) 

No 

Yes 

KDOV  Vis  LT  0.5  mi  (T+2) 

No 

Yes 

KDOV  Vis  LT  0.5  mi  (T+3) 

Yes 

Yes 

KDOV  Vis  LT  0.5  mi  (T+4) 

Yes 

Yes 

KDOV  Vis  LT  0.5  mi  (T+5) 

Yes 

Yes 

KDOV  Vis  LT  0.5  mi  (T+6) 

Yes 

Yes 

KWRI  Cig  LT  200  ft  (T+1) 

No 

Yes 

KWRI  Cig  LT  200  ft  (T+2) 

Yes 

Yes 

KWRI  Cig  LT  200  ft  (T+3) 

Yes 

Yes 

KWRI  Cig  LT  200  ft  (T+4) 

Yes 

Yes 

KWRI  Cig  LT  200  ft  (T+5) 

Yes 

Yes 

KWRI  Cig  LT  200  ft  (T+6) 

Yes 

Yes 

KWRI  Vis  LT  0.5  mi  (T+1) 

No 

No 

KWRI  Vis  LT  0.5  mi  (T+2) 

Yes 

Yes 

KWRI  Vis  LT  0.5  mi  (T+3) 

Yes 

Yes 

KWRI  Vis  LT  0.5  mi  (T+4) 

Yes 

Yes 

KWRI  Vis  LT  0.5  mi  (T+5) 

Yes 

Yes 

KWRI  Vis  LT  0.5  mi  (T+6) 

Yes 

Yes 

KADW  Cig  LT  200  ft  (T+1) 

No 

No 

KADW  Cig  LT  200  ft  (T+2) 

Yes 

No 

KADW  Cig  LT  200  ft  (T+3) 

No 

Yes 

KADW  Cig  LT  200  ft  (T+4) 

Yes 

Yes 

KADW  Cig  LT  200  ft  (T+5) 

Yes 

Yes 

KADW  Cig  LT  200  ft  (T+6) 

Yes 

Yes 

KADW  Vis  LT  0.5  mi  (T+1) 

No 

Yes 

KADW  Vis  LT  0.5  mi  (T+2) 

No 

Yes 

KADW  Vis  LT  0.5  mi  (T+3) 

No 

Yes 

KADW  Vis  LT  0.5  mi  (T+4) 

No 

Yes 

KADW  Vis  LT  0.5  mi  (T+5) 

Yes 

Yes 

KADW  Vis  LT  0.5  mi  (T+6) 

Yes 

Yes 

where  Y  is  the  dependant  or  response  variable,  in  this  case,  the  predietand.  Eaeh  x  is  a 
regressor  variable,  in  this  ease  the  predietors,  having  a  eorrelation  with  the  dependant 
variable.  Eaeh  ^^is  a  regression  eoeffieient,  caleulated  using  a  method  of  least  squares. 
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which  minimizes  the  error  between  the  predictand  (7)  and  eaeh  predietor  (x).  The 
response  variable  is  interpreted  as  the  expected  value  of  Y  for  any  speeifie  combination  of 
x’s,  the  regressor  variables  (Neter  et  al.  1990).  For  this  reason,  Y  is  interpreted  as  a 
probability  foreeast  for  future  oeeurrenee  or  non-oeeurrenee  of  a  partieular 
meteorologieal  event.  Finally,  £■  is  a  random  error  term. 

Predietors  are  added  to  the  regression  equations  based  on  a  mixed  stepwise 
regression  teehnique.  A  mixed  stepwise  teehnique  adds  regressor  variables  to  the 
predietive  equation  one  at  a  time  based  on  the  lowest  F-statistie.  As  more  variables  are 
added,  F-statisties  are  reeomputed  and  regressor  variables  may  be  removed  from  the 
equation  if  their  F-statistie  value  is  inereased  based  on  interaction  between  terms.  The 
major  advantage  to  using  a  mixed  stepwise  teehnique  over  a  forward  teehnique  is  the 
ability  to  remove  terms  whose  F-statistie  value  has  inereased  after  new  terms  are  added. 
This  takes  into  aeeount  interaetion  between  predietors,  thereby  minimizing  the  ehanees  of 
redundant  or  interaetive  predietors  remaining  in  the  final  predictive  equation. 

The  F-test  is  the  eriteria  for  adding  a  predietor  to  a  stepwise  regression 
(Montgomery  and  Runger  2003);  there  is  no  limit  to  the  number  of  predietors  in  eaeh 
equation,  as  long  as  they  meet  the  F-test  requirements  (whieh  are  set  stringently  at 
P<0.001).  The  F-statistie  is  ratio  between  the  sum  of  the  square  errors  of  the  eoefficients 
(P)  and  the  mean  square  error  (MSB)  of  the  dependant  variables  (Montgomery  and 
Runger  2003).  For  this  reason,  a  minimized  F-statistie  produees  a  model  with  less  error 
and  is  therefore  most  desirable.  Surfaee  and  upper  air  predietors  are  added  or  removed 
individually  to  eaeh  regression  equation  based  on  the  above  method. 
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Development  of  a  single  equation  for  one  predietand  at  a  speeifie  time  in  the  future 
oeeurs  as  follows  (Fig  5).  Predietor  sets  from  eaeh  network  surfaee  observing  station, 
upper  air  level,  and  elimatologieal  data  are  first  regressed  individually  against  the 
predietand  value  (e.g.  KDOV  CIG  LT  200  ft),  using  an  F-test  signifieanee  of  0.001.  Eaeh 
predietor  that  meets  the  striet  0.001  F-test  signifieanee  level  is  then  eombined  and 
regressed  against  the  predietand  onee  again,  using  a  mueh  looser  signifieanee  level  of 
F<0.1.  Sinee  all  terms  meet  the  original  F-statistie  eriteria,  the  larger  error  rate  is 
aeeepted  simply  to  remove  interaetive  or  redundant  terms.  This  proeedure  produees  one 
regression  equation  for  eaeh  predietand  at  eaeh  time  interval. 

To  produee  a  single  probabilistie  foreeast,  data  are  fed  into  the  regression  equations, 
whieh  prediet  a  value  between  0  and  1  that  ean  be  interpreted  as  a  foreeasted  probability 
of  oeeurrenee  (Visloeky  and  Fritseh  1997).  For  example,  the  oeeurrenee  of  a  eeiling  less 
than  200  ft  at  Dover  AFB  is  eoded  as  a  ‘  1’.  After  the  regression  equation  is  ealeulated, 
the  resulting  value  is  the  probability  of  the  200  ft  eeiling  existing  at  some  future  time.  A 
problem  with  linear  regression,  as  noted  by  Wilks  (1995)  is  that  in  rare  eases  probabilities 
ean  fall  outside  of  the  0  to  1  probability  range.  However,  this  is  not  a  signifieant  problem 
for  the  operational  meteorologist,  as  foreeast  values  outside  of  this  range  ean  be  rounded 
to  0  or  1  and  still  produee  aeeurate  foreeasts  (Wilks  1995).  All  statistieal  analysis  is 
performed  using  the  JMP  statistieal  software. 

After  predietive  equations  are  developed,  verifieation  is  required  for  eaeh  airfield 
for  eaeh  predietand  at  eaeh  time  period.  Verifieation  is  aeeomplished  on  an  independent 
data  set.  Eaeh  predietive  equation  is  eomputed  and  the  probability  foreeast  value  is 
eompared  with  the  eonditional  elimatology  foreeast  (onee  again  obtained  through 
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FIG  5.  Linear  Regression  Flowehart. 


AFCCC).  Foreeasts  are  produeed  for  a  one  year  independent  data  set  for  eaeh  fog 
oeeurrenee  using  both  the  observations-based  method  and  eonditional  elimatology.  This 
methodology  is  detailed  in  Figure  6. 
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Assessment  of  foreeast  aeeuraey  is  first  made  by  caleulating  the  mean  square  errors 
for  both  foreeast  methods.  The  MSB  is  the  averaged  squared  differenee  between  the 
foreeast  and  actual  event.  The  lower  the  MSB,  the  higher  the  accuracy  of  the  forecast.  A 
skill  score  was  then  calculated  per  Wilks  (1995)  to  show  the  percent  improvement  of  the 


BIG  6.  Verification  Blowchart. 
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observations-based  network  over  conditional  climatology. 

MSE  =  —  ■  (yfcgt  —  Yobs) 
n 

SS=l-t!!!2!2  (4) 

MSEcc 

Attempts  are  made  to  quantify  the  accuracy  of  the  forecast  in  an  operational 
environment.  The  first  category  evaluated  is  “bad  forecasts,”  which  is  defined  as  a 
forecast  with  an  absolute  error  greater  than  0.50— e.g.  more  than  50%  error  between  the 
forecast  and  observed  conditions.  A  second  assessment  is  made  of  “good  forecasts.” 
These  are  defined  as  having  an  absolute  error  of  less  than  0.30.  In  order  for  this 
methodology  to  be  implemented  operationally,  good  forecast  must  be  the  norm.  Forecast 
probabilities  between  30-60%  have  little  operational  value  to  the  operational  forecaster, 
regardless  of  the  MSB.  Mean  square  errors,  number  of  good  forecasts,  and  number  of 
bad  forecasts  are  computed  for  both  the  observations-based  network  and  conditional 
climatology  and  are  compared. 

This  chapter  details  the  methodology  employed  to  develop  a  forecast  tool  for 
accurate  prediction  of  fog  dissipation.  An  observations-based  network  consisting  of 
surface  observations,  upper  air  observations,  and  climatological  parameters  is  developed 
for  the  15  OWS.  A  multiple  linear  regression  model  is  implemented  for  each  of  the 
predictands  at  the  three  15  OWS  airfields. 

The  results  of  the  research  are  now  examined.  The  initial  goal  is  to  provide  a 
superior  forecast  technique  to  conditional  climatology.  To  accomplish  this,  36 
probabilistic  forecasts  are  developed  for  the  predictands  listed  in  Table  10.  The  next 
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chapter  details  the  suecesses  and  failures  of  the  observations-based  system  and  examines 
the  eharaeteristies  of  the  predictive  equations. 
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IV.  Analysis  and  Results 


4.1  Overview 

Predictive  equations  are  developed  using  the  aforementioned  statistical  technique 
for  two  predictands,  ceiling  less  than  200  ft  and  visibility  less  than  0.5  mi,  at  each 
airfield,  for  hourly  intervals  one  to  six  hours  in  the  future,  for  a  total  of  36  predictive 
equations.  Of  the  36  equations,  16  are  able  to  consistently  better  conditional  climatology 
when  analyzed  against  an  independent  data  set. 

This  section  provides  predictive  equations,  verification  statistics  on  the 
independent  data  set,  analysis  of  the  predictive  equations,  and  development  of  trends. 
Analysis  on  the  successes  and  limitations  of  the  technique  are  examined.  Section  4.2 
details  the  results  obtained  at  Dover  AFB,  with  the  following  sections  looking  at 
McGuire  AFB  and  Andrews  AFB  individually.  After  each  individual  airfield  is 
examined,  underlying  problems  persistent  throughout  the  forecast  technique  are 
thoroughly  examined.  Finally,  examples  of  application  of  the  methodology  to  real  world 
situations  are  examined. 

4.2  Dover  AFB 

4.2.1  Ceiling  Less  than  200 ft.  Mean  square  errors  are  calculated  using  both  conditional 
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climatology  and  the  observations-based  network  for  an  independent  one  year  sample  of 
data  consisting  of  approximately  45  oceurrenees  and  are  graphed  in  Figure  7a.  Note  that 
the  MSEs  for  the  observations-based  system  are  superior  to  conditional  climatology  at 
each  forecast  time.  Note  also  deerease  in  the  MSEs  with  time.  While  this  is  eounter 
intuitive,  the  set  up  of  the  experiment  results  in  this  trend.  Oecurrences  of  eeilings  below 
200  ft  and  visibility  less  than  0.5  are  rare  and  typically  short-lived  events,  therefore,  there 
is  less  variability  in  the  observed  conditions  at  the  later  hours  (e.g.  the  observed  condition 
is  coded  as  a  “0”  representing  non-oecurrence  for  a  majority  of  the  events).  The 
observations-based  network  and  conditional  climatology  typically  produce  low  foreeast 
probabilities  at  these  hours  resulting  in  the  decreased  MSE  with  time.  Three  measures 
of  foreeast  aceuracy,  as  discussed  in  the  methodology  seetion,  are  also  shown  in  Eigures 
7b-d. 

Eirst,  skill  scores  for  each  forecast  time  are  illustrated,  with  positive  values 
indicating  an  improvement  over  elimatology  of  the  observations-based  network.  All  six 
foreeast  times  show  significant  improvement  over  conditional  climatology.  Eorecast 
improvements  range  from  1.3%  to  50%,  with  an  average  improvement  of  15.8%.  The 
bad  foreeast  category  (Eig.  7d)  shows  that  the  observations-based  network  is  equal  to  or 
superior  to  eonditional  climatology  at  all  six  hourly  forecasts.  Bad  forecasts  range  from 
2-18  per  hour  for  the  observations-based  network  and  2-26  per  hour  for  eonditional 
elimatology.  The  extreme  case  demonstrates  conditional  climatology  can  be  inaceurate 
on  over  50%  of  the  foreeasts.  The  final  foreeast  measure,  good  foreeasts  (Eig.  7c),  again 
demonstrates  a  significant  improvement  over  conditional  climatology.  The  only  ease 
where  conditional  climatology  is  better  is  at  the  6-hour  point  where  it  is  favored  due  to 
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the  fact  that  200  ft  ceilings  are  typically  a  short-lived  event.  On  average,  the 
observations-based  network  produces  15%  more  good  forecasts  at  each  hourly  interval. 

Predictive  equations  for  Dover  AFB  ceilings  are  available  in  Table  11.  The 
number  of  predictors  selected  for  each  equation  varies  from  10  at  the  earliest  times  to  16 
at  the  latest  forecast  hours.  Standardized  beta  (std  beta)  values  are  included  in  order  to 
gauge  the  relative  strength  of  each  predictor  on  each  model.  The  highest  absolute  value 
among  the  standardized  betas  have  the  most  influence  on  the  forecast.  Predictors  are  now 
examined  for  trends. 


MSE  for  KDOV  Ceiling  Less  Than  200  ft 
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FIG.  7.  Forecast  Statistics  for  KDOV  Ceiling  Less  than  200  ft.  Panel  (a)  displays  MSE, 
panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad  forecasts 
for  each  hour. 
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TABLE  11.  Forecast  Equations  for  KDOV  Ceiling  Eess  than  200  ft.  The  first  3  hours  of  the  forecast  (a)  and  the  last  three  (b). 
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500MBRH>50%  -0.015  -0.133  KEZF  VIS  <  4800  M  -0.115  -0.120  KFRG  CIG  HI  <1000  FT  0.101  0.133 

0.167  0.119  KCJR  SKY  (FEW/SCT)  0.116  0.090  KLFI  CIG  HT  <200  FT  0.108  0.088 

KFRG  SKY  (FEW/SCT)  -0.214  -0.087  KPTB  SKY  (OVC)  -0.123  -0.177 

b)  700  MB  T  -0.014  -0.138  KPTB  CIG  HT  <200  FT  0.097  0.085 

500MBRH>50%  0.198  0.156  500  MB  T  -0.010  -0.110 


A  successful  trend  is  the  KDOV  SKY  (OBSCURED)  predietor  showing  up  in  the 
first  few  hours  of  the  foreeast.  This  is  important  beeause  it  is  indieates  a  foggier  situation 
that  may  persist  longer  into  time.  The  standardized  beta  values  deerease  with  time, 
indieating  that  this  eondition  has  less  impaet  on  the  foreeast  as  time  advanees.  The  least 
influential  region  for  predietors  eomes  from  stations  due  west  of  Dover  AFB.  The  early 
time  periods  demonstrate  a  uniformity  of  direetion  among  the  seleeted  network  stations, 
however,  as  time  advanees,  stations  to  the  southwest  beeome  more  dominant  in  the 
equations.  These  trends  are  expected,  sinee  weather  typically  moves  from  west  to  east  in 
mid  latitudes.  All  geographieal  direetions  infiuenee  the  earlier  hours  as  the  situation  is 
stagnant  under  typieally  stable,  light  wind  environment. 

Another  sueeessful  trend  is  the  presenee  of  at  least  one  easterly  wind  predietor 
(bold.  Table  11),  which  indicates  a  higher  likelihood  of  eeilings  remaining,  due  to  low- 
level  moisture  adveetion  from  the  Atlantic  Ocean.  These  predietors  maintain  their 
strength  throughout  the  six  foreeast  hours.  The  faet  that  the  only  upper  air  predictors  are 
temperature  and  moisture  predietors,  not  wind  predietors,  indieates  the  model  deteets 
eeilings  dissipated  through  lifting  proeesses  better  than  advective  proeesses. 

Of  eoneern  is  the  large  number  of  oeeurrenees  of  unrestrieted  visibility  as  a 
dissipation  predietor.  While  unrestricted  visibilities  reported  at  stations  elose  to  Dover 
AFB  indieate  dissipation,  stations  further  away  indieate  the  persistenee  of  the  eeiling. 
While  visibility  and  eeiling  eertainly  have  a  eorrelation,  it  is  a  eoneern  that  this  is  one  of 
the  most  influential  predietors,  espeeially  in  the  early  foreeasts. 

Overall,  the  eeiling  foreeasts  for  Dover  AFB  are  very  sueeessful  in  defeating 
eonditional  elimatology,  produeing  low  MSEs,  a  large  number  of  good  foreeasts  and  a 
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small  number  of  bad  forecasts.  Section  4.2.2  looks  at  the  visibility  less  than  0.5  mi 
forecast  for  Dover  AFB. 

4.2.2  Visibility  Less  than  0.5  mi.  The  forecasts  for  visibility  at  Dover  AFB  are  not  as 
successful  as  the  ceiling  forecasts.  Mean  square  errors  are  calculated  for  each  forecast 
method  and  are  displayed  in  Figure  8a.  The  MSEs  for  the  observations-based  network, 
while  still  relatively  low,  are  higher  than  those  for  conditional  climatology  for  each 
forecast  time.  While  the  observations  network  may  still  produce  accurate  forecasts,  in 
this  case  conditional  climatology  is  more  accurate,  when  verified  against  the  independent 
data  set. 

The  skill  scores  and  forecast  statistics  determine  if  there  are  any  advantages  to  the 
observations-based  network.  Figure  8b  shows  an  average  39%  decrease  in  skill  for  the 
observations-based  network.  The  observations-based  network  averages  8  more  bad 
forecasts  (out  of  a  sample  size  of  approximately  150)  at  each  hour.  While  not  a  positive 
result,  this  is  a  better  percentage  of  bad  forecasts  than  conditional  climatology  exhibits  on 
the  ceiling  forecasts.  A  positive  sign  for  the  observations-based  network  lies  in  the  good 
forecasts.  In  the  first  hour  forecast,  the  observations-based  network  produces  19  more 
good  forecasts,  and  in  the  second  hour  they  are  even.  Of  the  100  times  fog  persisted  into 
the  one  hour  forecast,  the  observations-based  network  produced  a  good  forecast  78  times. 
This  is  positive  because  the  first  couple  of  hours  are  typically  the  most  critical  and  often 
highly  variable. 

Upon  closer  examination,  the  bad  forecasts  (e.g.  higher  MSEs)  can  be  explained 
by  the  data.  The  first  problem  is  overestimating  non-occurrences  of  fog  through  all  six 
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MSE  for  KDOV  Visibility  Less  Than  0.5  mi 


Skill  Scores  for  KDOV  Visibility  Less  Than  0.5  mi 


FIG.  8.  Forecast  Statistics  for  KDOV  Visibility  Less  than  0.5  mi.  Panel  (a)  displays 
MSE,  panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad 
forecasts  for  each  hour. 


hours  of  the  forecast.  This  is  illustrated  by  the  fact  that  the  average  predicted  value  of 
occurrence  at  the  one -hour  point  is  76%.  This  results  in  the  high  number  of  misses  on  the 
situations  where  fog  is  not  present,  essentially  driving  up  the  false  alarm  rate  and  MSE. 
This  problem  of  overestimating  the  probability  of  occurrence  is  magnified  in  the  later 
forecasts  as  the  actual  chance  of  fog  is  more  likely  to  be  zero.  The  second  reason  for  the 
higher  MSE  is  a  probability  forecast  of  1  or  0  being  produced  and  the  opposite  value 
being  observed.  This  problem  is  magnified  by  the  set  up  of  the  experiment.  That  is  if  the 
observations  network  is  not  accurately  forecasting  fog  for  a  particular  event,  it  is  very 
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likely  to  miss  the  oeeurrenee  of  fog  at  all  time  intervals.  Therefore  the  MSB  is  inereased 
for  eaeh  foreeast  hour.  Examination  of  the  predietive  equations  as  shown  in  Table  12, 
eould  offer  a  possible  physieal  explanation  for  the  inferior  foreeasts. 

Analyses  of  the  equations  reveal  two  signifieant  faetors  most  strongly  influeneing 
a  dissipation  foreeast.  This  first  is  in  the  short-term  (one  to  two  hours  and  to  a  lesser 
extent  the  third  hour  foreeast)  and  oeeurs  at  the  Georgetown,  DE  observing  station 
(KGED,  bold  in  Table  12).  Two  to  three  values  per  foreeast  prediet  dissipation  among 
these  predictors  (the  negative  estimates  in  the  equations).  Of  particular  note  is  the 
influence  of  broken  sky  conditions  at  KGED  having  one  of  the  strongest  correlations  with 
dissipation  (based  on  standardized  beta  values)  at  each  of  the  first  two  hours. 

Examination  of  the  independent  data  set  reveals  that  on  occasions  with  inaccurate 
dissipation  forecasts,  these  conditions  are  not  present  at  KGED,  which  in  turn  keep  the 
probability  forecast  at  a  higher  level  than  it  should  be.  With  dissipation  probability  so 
closely  tied  to  this  one  station,  especially  in  the  early  forecasts,  slight  changes  in 
meteorological  parameters  (e.g.  sky  condition  broken  versus  overcast)  have  significant 
impact  on  producing  an  inaccurate  dissipation  forecast. 

The  second  predictor  subset  leading  to  a  dissipation  forecast  is  wind  direction 
between  south  and  northwest  (italicized  in  Table  12).  Elnfortunately,  this  condition  is  not 
present  in  a  majority  of  cases  where  dissipation  occurs.  This  is  found  throughout  all 
forecasts  and  is  related  to  these  particular  wind  directions,  regardless  of  station.  This 
finding  implies  that  this  method  has  difficulty  detecting  dissipation  unless  the  correct 
wind  sector  predictors  are  present. 
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TABLE  12.  Forecast  Equations  for  KDOV  Visibility  Eess  than  0.5  mi.  The  first  3  hours  of  the  forecast  (a)  and  last  three  (b). 
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Finally,  the  equations  developed  for  the  first  two  forecast  hours  exhibit  underlying 
problems  in  the  model  data  set.  Note  from  Table  10  these  two  data  sets  have 
distributions  of  errors  which  are  not  normal.  This  statistical  problem  helps  explain  some 
of  the  inaccuracies  the  first  two  models  demonstrate. 

Overall  the  forecasts  for  Dover  AFB  show  promise,  especially  in  the  ceiling 
category.  Although  conditional  climatology  has  a  lower  MSB,  there  are  advantages  to 
the  visibility  forecasts,  especially  when  the  good  forecasts  in  the  very  short  term  are 
considered.  The  predictive  equations  for  McGuire  AFB  are  examined  next. 

4. 3  McGuire  AFB 

4.3.1  Ceiling  Less  than  200 ft.  Mean  square  errors  for  the  six  forecast  hours  are 
illustrated  in  Figure  9a.  The  MSEs  are  smaller  for  the  observations-based  network  for  all 
six  forecast  hours.  Examination  of  the  skill  scores  (Eig.  9b)  show  improvement  over 
conditional  climatology  for  each  of  the  six  forecast  hours  ranging  from  3%-49%  with  an 
average  improvement  of  22.5%.  The  number  of  good  forecasts  (Eig.  9c)  is  consistent 
with  both  forecast  methods,  averaging  about  32  good  forecasts  out  of  approximately  55 
samples.  Aside  from  the  skill  score,  the  other  distinctive  advantage  of  the  method  is  the 
minimization  of  bad  forecasts  (Eig  9d).  Conditional  climatology  produces  an  average  6% 
more  bad  forecasts  per  hour. 

An  examination  of  the  forecast  equations  (Table  13)  leads  one  to  conclude  which 
variables  produce  the  most  accurate  forecasts.  Eirst,  as  the  forecast  advances  into  the 
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FIG.  9.  Forecast  Statistics  for  KWRI  Ceiling  Less  than  200  ft.  Panel  (a)  displays  MSE, 
panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad  forecasts 
for  each  hour. 


future,  the  climatological  factors  are  among  the  stronger  predictors  influencing  the 
equations.  This  is  expected,  since  in  general,  as  the  forecast  time  increases,  the  relevance 
of  the  current  observation  decreases.  Second,  the  geographic  region  of  the  predictor 
station  relative  to  McGuire  AFB  has  a  significant  influence,  especially  as  time  increases. 
Stations  from  the  southwest  and  west  dominate  the  predictive  equations,  especially  in  the 
later  stages.  This  is  consistent  with  the  findings  for  the  ceiling  at  Dover  AFB,  where  the 
southwest  and  west  stations  dominate  the  predictive  equations  at  the  longer  time  periods. 
Of  the  96  predictors  in  the  six  equations,  wind  predictors  dominate  as  43  (indicated  in 
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TABLE  13.  Forecast  Equations  for  KWRI  Ceiling  Eess  than  200  ft.  The  first  3  hours  of  the  forecast  (a)  and  the  last  three  (b). 
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bold  in  Table  13)  are  based  on  wind  speed  or  direetion.  Wind  speed  and  dew  point 
depression,  have  a  signilieant  effeet  on  the  predietion  of  dissipation,  espeeially  in  the 
earlier  foreeast  hours.  These  are  positive  signs,  sinee  they  are  two  of  the  most  important 
meteorologieal  parameters  indieating  dissipation.  Also  of  note  is  the  higher  number  of 
upper  level  predietors  as  eompared  to  the  Dover  AFB  eeiling  foreeast  equations. 

A  final  faetor  indieating  the  strength  and  aeeuraey  of  this  model  is  the  laek  of 
dependenee  on  erroneous  predietors.  Whereas  the  Dover  AFB  eeiling  foreeast  appears 
overly  dependant  on  visibility  predietors,  relatively  few  are  ineluded  in  MeGuire  AFB’s 
eeiling  foreeast.  The  equations  are  heavily  infiueneed  by  wind,  moisture,  eeiling,  and  sky 
eondition  predietors  as  would  be  expeeted  with  a  eeiling  foreeast  model.  Next  examined 
is  the  MeGuire  AFB  visibility  foreeast  that  exhibits  some  advantages  despite  not  being  as 
strong  as  the  eeiling  foreeast. 

4.3.2  Visibility  Less  than  0.5  mi.  Mean  square  errors  for  both  foreeast  teehniques  are 
shown  in  Figure  10a.  Half  of  the  foreeasts  of  the  observations-based  network  are 
superior,  while  the  other  three  favor  eonditional  elimatology.  The  three-hour  foreeast 
shows  the  largest  differenee  in  MSB,  likely  due  to  the  large  differenee  in  number  of  bad 
foreeasts  between  the  two  teehniques  (Fig.  lOd).  Removing  the  three-hour  foreeast 
leaves  the  overall  average  MSEs  almost  even.  Figure  10b  shows  a  wide  variety  of  range 
among  the  skill  seores,  from  an  improvement  of  34%  in  the  observations-based  network 
to  a  deerease  in  skill  of  14%,  with  an  average  inerease  of  3%.  The  observations-based 
network  shows  a  higher  number  of  good  foreeasts  for  the  first  half  of  the  foreeast  period 
while  eonditional  elimatology  beeomes  superior  at  the  later  foreeast  hours  (Fig.  lOe). 
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FIG.  10.  Forecast  Statistics  for  KWRI  Visibility  Less  than  0.5  mi.  Panel  (a)  displays 
MSE,  panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad 
forecasts  for  each  hour. 


The  equations  for  each  forecast  are  in  Table  14.  There  are  two  significant 
differences  when  comparing  the  one-hour  forecast  equation  to  the  2,  3,  and  4  hour 
forecast  equations,  which  are  the  successful  equations.  The  first,  and  least  significant,  is 
the  absence  of  any  climatological  parameters.  Climatological  parameters  are  not  believed 
to  be  the  main  source  of  error  in  equation  one  because  they  are  typically  more  influential 
on  the  later  forecasts.  The  second  major  ingredient  missing  is  the  presence  of  any  dew 
point  or  dew  point  depression  parameters  from  the  one-hour  equation.  At  least  one  or 
more  of  these  predictors  generally  appear  in  all  of  the  other  equations  and  are  the  most 
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important  factors  in  fog  formation  or  dissipation.  This  is  believed  to  be  the  reason  for  the 
low  skill  seore  (Fig.  10b)  in  the  one-hour  foreeast.  However,  there  are  still  advantages  to 
the  one  hour  foreeast,  ineluding  its  produeing  81  good  foreeasts  (Fig.  lOe)  out  of  150, 
signifieantly  more  than  the  66  issued  by  eonditional  elimatology. 

The  five  and  six-hour  equations  also  have  negative  skill  seore,  however  moisture 
and  elimatologieal  parameters  are  ineluded  in  them.  The  reason  for  the  deerease  in  skill 
seore  in  this  ease  is  over  predietion  of  the  event.  While  bad  foreeasts  (Absolute  Error  > 
50%)  are  eomparable  with  eonditional  elimatology,  foreeast  values  between  30%-50% 
produee  a  large  foreeast  error.  Due  to  eonditional  elimatology  generally  eonverging  to 
zero  at  the  later  hours,  it  produees  lower  MSEs  in  the  five  to  six  hour  period.  Eor 
example,  in  the  five  hour  foreeast,  the  average  predieted  value  using  the  observations- 
based  network  is  21%  while  the  eonditional  elimatology  average  is  17%.  The  lower 
average  foreeast  value  produees  the  lower  MSE  when  over  75%  of  the  observations  are 
for  non-oeeurrenee  of  the  event. 

The  most  influential  geographieal  region,  as  evideneed  in  the  2,  3,  and  4-hour 
foreeasts,  is  southwest.  This  is  another  faetor  in  the  degradation  of  the  foreeast  skill  in 
the  5  and  6-hour  foreeasts.  A  more  omni-direetional  sampling  of  station  data  is  evident  at 
the  longer  hours,  whieh  is  eontrary  to  sueeessful  foreeasts  whieh  are  typieally  heavily 
based  on  the  westerly  and  southwesterly  observations.  The  laek  of  influenee  of  the 
dominant  flow  on  the  later  hour  foreeasts  likely  results  in  the  deereased  skill  seores  at 
these  times. 

Note  that  the  upper  air  predietors  are  some  of  the  weakest  throughout  the  foreeast 
period,  while  the  elimatologieal  predietors,  whieh  appear  in  five  of  the  six  equations,  are 
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TABLE  14.  Forecast  Equations  for  KWRI  Visibility  Eess  than  0.5  mi.  The  first  3  hours  of  the  forecast  (a)  and  last  three  (b). 
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some  of  the  strongest.  Other  than  the  climatologieal  predietors,  there  are  no  distinct 
trends  among  the  meteorological  predictor  selection. 

Of  the  12  predictive  equations  developed  for  McGuire  AFB,  nine  show  superior 
forecast  skill  to  conditional  climatology.  The  unsuccessful  predictive  equations  produce 
good  forecasts,  with  slight  biases.  Overall,  the  application  of  the  surface-based  observing 
method  for  McGuire  AFB  is  considered  a  success.  Next,  evaluation  of  the  forecast 
equations  for  Andrews  AFB  is  examined. 

4. 4  Andrews  AFB 

4.4.1  Ceiling  Less  than  200 ft.  The  graph  of  MSEs  verified  on  the  independent  sample  of 
data  is  shown  in  Figure  11a.  It  is  evident  in  all  Andrews  AFB  cases  that  the  MSE  for 
conditional  climatology  is  smaller  than  that  of  the  surface-based  observations  system. 
Although  the  one  and  six-hour  forecasts  are  competitive,  the  other  four  clearly  are  not. 
The  main  reason  for  the  variability  in  MSEs  is  a  small  sample  size.  While  the  Dover 
AEB  and  McGuire  AEB  forecast  parameters  have  approximately  50  observations  for 
verification,  Andrews  only  has  approximately  30.  There  are  a  couple  of  reasons  for  this. 
Eirst,  the  number  of  occurrences  of  ceiling  below  200  ft  at  Andrews  AEB  is  significantly 
lower  than  the  other  two  stations.  Dover  AEB  and  McGuire  AEB  have  approximately 
400  observations  of  fog,  while  Andrews  AEB  has  only  about  200.  The  second  reason  is 
the  non-reporting  of  ceilings.  As  previously  mentioned,  forecasts  are  not  verified  if 
ceiling  data  is  not  available  at  the  forecast  hours.  This  is  especially  a  problem  at 
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MSE  for  KADW  Ceiling  Less  Than  200  ft 


Skill  Scores  for  KADW  Ceiling  Less  Than  200  ft 


FIG  11.  Forecast  Statistics  for  KADW  Ceiling  Less  than  200  ft.  Panel  (a)  displays 
MSE,  panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad 
forecasts  for  each  hour. 


Andrews  AFB,  where  the  occurrence  of  reduced  visibility  often  resulted  in  no  sky 
condition  being  reported. 

A  second  problem  in  the  Andrews  AFB  data  set  is  the  large  amount  of  missing 
upper  air  data.  Section  3.2.3  shows  that  12%  of  the  upper  air  data  is  missing  for  Andrews 
AFB,  much  higher  than  the  other  two  data  sets.  This  results  in  fewer  upper  air  predictors 
being  selected  for  each  forecast.  A  final  problem  to  note  in  the  Andrews  AFB  data  set  is 
the  failure  to  meet  the  necessary  assumptions  required  for  linear  regression.  Table  10  in 
Chapter  III  illustrates  that  the  first  three  forecast  hours  exhibit  non-normality. 
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heteroscedasticity,  or  both.  The  linear  regression  model  is  probably  not  applieable  to  this 
data  set  and  has  a  strong  influenoe  on  the  exeessively  high  MSEs  in  the  first  three  hours. 
The  average  skill  seore  for  the  six  foreeasts  is  -46%  (Fig.  1  lb).  The  one-hour  foreeast  is 
most  eompetitive  with  eonditional  elimatology,  it  has  a  skill  seore  of -1%,  but  it  offers 
fewer  bad  forecasts  than  conditional  climatology  and  almost  the  same  MSE  (Fig.  1  la  & 
lid).  Aside  of  the  inherent  problems  with  the  data,  an  analysis  of  the  equations  produced 
by  the  observations-based  network  (Table  15)  provides  some  insight  into  why  the  results 
from  Andrews  AFB  base  differ  from  Dover  AFB  and  McGuire  AFB. 

An  analysis  of  the  equations  reveals  significant  findings.  First,  the  number  of 
predictors  varies  from  4-15  throughout  the  six  equations.  This  is  much  lower  than  the 
previous  equations,  which  typically  vary  from  about  10-25  predictors.  The  small  data  set 
does  not  adequately  detect  the  dominant  atmospheric  conditions  resulting  in  dissipation, 
therefore,  the  regression  is  unable  to  select  ample  predictors.  The  smaller  number  of 
predictors  (which  as  stated  before  are  mostly  binary),  unfortunately  are  not  representative 
enough  of  the  atmosphere  surrounding  the  station.  Variations  in  the  atmospheric 
conditions  that  could  lead  to  accurate  dissipation  predictions  are  therefore  not  well 
represented  by  the  model. 

Upper  air  and  climatological  predictors  are  also  omitted  from  the  equations. 

These  are  often  the  driving  forces  in  the  previous  successful  forecast  models.  Uniformity 
of  station  selection  is  also  evident  in  the  later  forecasts.  Successful  models  typically  have 
a  majority  of  the  stations  used  in  the  prediction  from  the  southwest  and  west  in  the  later 
periods.  This  model  has  an  omni-directional,  or  uniform  geographic  distribution  in  the 
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TABLE  15.  Forecast  Equations  for  KADW  Ceiling  Eess  than  200  ft.  The  first  3  hours  of  the  forecast  (a)  and  the  last  three  (b). 
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KCHO  DEWPOINT  DEPRESSION  0.065  0.177  KWWD  WIND  DIR  (<90)  0.180  0.218  KTTN  VIS  <  9999  M  0.135  0.086 

KPNEPRECIP  -0.253  -0.219  SEC  WIND  N  0.125  0.160  KBLM  CIG  HT  <1000  FT  0.I7I  0.182 

KPNE  WIND  SPEED  0.023  0.226  1 000  MB  WIND  N  0.I5I  0.184 


later  forecast  hours,  which  is  also  troublesome  in  the  Dover  AFB  and  McGuire  AFB 


cases  where  it  appears.  The  predictors  selected  by  the  model  do  not  appear  ambiguous. 
That  is,  the  model  focuses  mainly  on  sky  conditions  and  ceiling  predictors  and  veers 
away  from  the  visibility  predictors.  Overall,  it  is  believed  that  with  a  larger,  normally 
distributed  sample  and  some  fine  tuning  this  model  could  be  made  to  accurately  represent 
the  forecast  conditions.  Examined  next  is  the  visibility  forecast  for  Andrews  AFB. 

4.4.2  Visibility  Less  than  0.5  mi.  The  MSE  (Eig.  12a)  for  the  one-hour  forecast  is 
superior  in  the  observations-based  system;  however,  it  is  the  only  case  that  illustrates  the 
superiority  of  this  methodology.  The  MSEs  are  significantly  better  for  conditional 
climatology  for  the  other  five  forecast  hours.  This  is  also  evident  in  the  skill  scores  (Eig. 
12b)  where  the  average  decrease  in  skill  score  of  the  observations-based  network  is  45%. 
The  other  forecast  statistics  (Eig  12c,  12d)  significantly  favor  conditional  climatology, 
with  the  observations  network  competitive  on  only  a  few  occasions. 

Once  again,  it  is  believed  that  the  largest  problem  for  the  Andrews  AEB  visibility 
forecasts  is  the  data  itself  The  number  of  occurrences  is  about  one  third  the  number  of 
occurrences  for  both  Dover  AEB  and  McGuire  AEB.  Also  note  from  Chapter  III  (Table 
7)  that  the  first  four  forecast  hours  exhibit  errors  that  are  not  normally  distributed.  This 
violates  the  assumptions  required  for  linear  regression.  This  is  the  main  reason  for  the 
inadequacies  in  the  Andrews  AFB  forecast. 

Examining  the  equations  (Table  16)  provides  further  insight  on  the  problems  with 
the  visibility  forecast.  The  best  forecast  time,  at  one  hour,  has  relatively  few  predictors, 
but  is  geographically  uniform  around  the  station.  The  model  remains  consistent  in  its 
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Good  Forecasts  for  KADW  Visibility  Less  than  0.5  Bad  Forecasts  for  KADW  Visibility  Less  than  0.5  mi 

mi 


FIG  12.  Forecast  Statistics  for  KADW  Visibility  Less  than  0.5  mi.  Panel  (a)  displays 
MSB,  panel  (b)  shows  skill  scores,  panel  (c)  good  forecasts,  and  panel  (d)  shows  bad 
forecasts  for  each  hour. 


selection  of  parameters  throughout  the  six  hours,  and  also  appears  to  hone  in  on  the 
meteorological  processes  that  lead  to  dissipation.  For  example,  at  least  one  moisture 
parameter  is  present  in  all  six  forecasts  (bold  in  Table  16)  and  the  standard  beta  values 
show  that  they  are  one  of  the  most  influential  factors  predicting  dissipation.  Another 
noticeable  trend  is  the  strong  presence  of  northeasterly  to  southeasterly  winds,  which 
would  intuitively  indicate  moisture  advection  from  the  Atlantic  Ocean.  In  the  first  hour’s 
forecast,  the  easterly  wind  term  indicates  dissipation,  while  the  other  five  forecasts 
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TABLE  16.  Forecast  Equations  for  KADW  Visibility  Eess  than  0.5  mi.  The  first  3  hours  of  the  forecast  (a)  and  last  three  (b). 
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remain  consistent  in  indicating  the  presence  of  reduced  visibility  with  these  wind 
conditions.  This  discrepancy  is  left  unexplained  in  the  statistical  analysis. 

The  later  forecast  hours  again  show  neither  climatological  predictors,  nor  a 
dependence  on  westerly  or  southwesterly  stations  for  predictors.  These  could  be  two 
contributing  factors  to  the  lack  of  skill  shown  later  in  the  period. 

Although  the  Andrews  AFB  forecasts  overall  fail  to  verify  well  on  independent 
data,  there  are  positive  trends,  such  as  the  lack  of  erroneous  predictor  selection  indicating 
that  success  could  be  achieved  in  the  future.  The  small  number  of  predictors,  due  to  the 
smaller  data  set,  has  a  large  influence  on  the  negative  skill  scores  of  the  model. 

However,  the  most  significant  problem  in  the  models  is  the  failure  of  the  data  set  to  meet 
the  assumptions  required  for  linear  regression.  Underlying  problems  throughout  the 
models  are  examined  next. 

4.5  Underlying  Model  Problems.  There  are  inherent  problems  that  persist  through  all  the 
models  that  are  addressed.  The  first  is  the  linearity  of  the  model.  In  linear  regression,  an 
R  value  is  calculated,  this  is  the  coefficient  of  determination  and  represents  the  adequacy 
of  fit  of  the  regression  model,  or  the  variability  in  the  data  accounted  for  by  the  model 
(Montgomery  and  Runger  2001).  The  R  values  range  from  0  to  1  with  1  being  a  perfect 
fit.  The  R  values  for  the  forecast  models  presented  in  this  work  typically  are  between 
0.25  and  0.6.  Although  these  values  may  seem  low,  they  are  not  of  major  concern  for  the 
following  reasons.  First,  R  values  can  be  erroneously  inflated  by  adding  an  excessive 
amount  of  terms  (Montgomery  and  Runger  2001).  This  study  keeps  with  the  strict  F-test 
criteria  to  minimize  the  number  of  predictors  added  to  the  equation.  Second,  this  is  a 
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linear  model  of  the  atmosphere,  an  entity  that  usually  behaves  non-linearly,  making  high 
R  values  diffieult  to  obtain.  Third,  aeeording  to  Montgomery  and  Runger  (2003)  large 
R  values  do  not  neeessarily  imply  aeeurate  predietions  of  future  events.  Finally,  as 
Wilks  (1995)  points  out  the  goal  of  the  statistieal  methodology  is  to  produee  an  aeeurate 
foreeast,  not  a  perfeet  model.  Therefore,  this  work  foeuses  on  developing  aeeurate 
dissipation  foreeasts,  regardless  of  the  model’s  linear  relationship  to  the  data  set. 

A  seeond  issue  is  the  absenee  of  sky  eonditions  in  observations.  Although 
replaeement  data  are  used  to  fill  in  the  missing  values,  sky  eondition  is  never  used  as  a 
predietand  value.  This  reduees  the  sample  size  eonsiderably  for  both  equation 
development  and  verifieation,  espeeially  in  the  Andrews  AFB  data  set. 

Third,  this  system  is  designed  independent  of  time  of  day  so  that  one  foreeast 
equation  eould  suffiee  for  the  airfield.  This  varies  from  Visloeky  and  Fritseh’s  (1997) 
work  that  foeused  on  two  short  time  intervals  for  their  foreeasts.  Breaking  the  data  into 
intervals  (2-4  hours)  would  be  more  representative  of  diurnal  differenees  in  weather 
patterns  and  eould  produee  improved  equations  for  fog  dissipation.  Furthermore,  it 
eliminates  some  of  the  advantages  eonditional  elimatology  has.  That  is,  as  fog  is 
primarily  a  noeturnal/early  morning  event  eonditional  elimatology  has  an  advantage  in 
the  afternoon  oases  in  that  probabilities  deorease  to  zero  muoh  faster,  typioally  mirroring 
the  observations.  Many  of  the  over  foreoasting  ooourrenoes  ooour  in  late  morning  or 
afternoon  situations,  where  the  observations-based  network  is  working  regardless  of  time. 
While  the  observations  network  is  based  on  the  physioal  meohanisms  that  lift  the  fog, 
different  prooesses  are  at  work  in  the  overnight  versus  daylight  hours.  Breaking  into  time 
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blocks  could  allow  the  equations  to  piek  different  predietors,  e.g.  different  physieal 
meehanisms  for  afternoon  versus  noetumal  events. 

A  final  remark  is  the  impaet  automated  observing  systems  employed  by  the  NWS 
had  on  this  study.  Visloeky  and  Fritseh’s  (1997)  original  work  used  abundant  full  time 
manual  reporting  stations  along  the  east  eoast,  many  of  whieh  no  longer  exist.  The 
disadvantage  of  the  ASOS  observations  is  missing  data  due  to  staggered  implementation 
and  spin-up  times.  These  two  faetors  eause  many  of  the  observations  from  ASOS 
stations  to  be  strongly  infiueneed  by  missing  data  teehniques.  Examples  of  the  how  this 
teehnique  is  applied  to  real  world  events  follow. 

4.6  Applications  of  the  Observations-Based  Method.  To  implement  this  method 
operationally,  eaeh  equation  is  ealeulated  to  produee  a  probability  foreeast.  Foreeast 
probabilities  range  from  ‘1’  representing  a  100%  ehanee  of  oeeurrenee  of  the  event  (e.g. 
eeiling  below  200  ft),  to  ‘0’  representing  a  0%  ehanee  of  the  event  oeeurring. 

Table  17  provides  an  applieation  of  this  proeess  to  produee  a  predietion  foreeast 
for  1-hour  in  the  future  at  Dover  AFB.  The  first  eolumn  represents  the  predietor  labels 
(the  physieal  meaning  of  the  response  variable)  with  the  regression  eoeffieients  (fi)  in  the 
seeond  eolumn.  The  third  eolumn  is  the  aetual  or  observed  value  of  the  predietor  (x). 
Notiee  that  a  majority  of  the  predietors  are  binary  values,  eoded  as  either  a  ‘0’  for  non- 
oeeurrenee  or  a  ‘  1’  for  oeeurrenee.  The  values  are  multiplied,  and  then  summed,  per  Eq. 
2  (ineluding  the  intereept,  fo)  to  produee  the  foreeast  of  76.9%  probability  of  oeeurrenee. 
Probabilities  are  ealeulated  for  eaeh  of  the  six  hourly  foreeasts  for  eaeh  of  the  two 
network  versus  eonditional  elimatology  on  aetual  observed  events.  The  first  event  took 
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Table  17.  Example  Predictive  Equation. 


Predictor 

Estimate 

(P) 

Observed 

Value 

Equation  Result 

Intercept: 

0.6999357 

0.6999357 

700  MB  T 

-0.018243 

4.6 

-0.0839178 

KWWD  CIG  HT  <3000  FT 

-0.4099 

0 

0 

KWWD  VIS  >=  9999  M 

-0.11243 

0 

0 

KSBY  VIS  >=  9999  M 

-0.120956 

0 

0 

KILG  VIS  >=  9999  M 

-0.11369 

0 

0 

KILG  WIND  DIR  (<235) 

-0.494969 

0 

0 

KMIV  WIND  DIR  (<90) 

0.2029688 

0 

0 

KGED  CIG  HT  <3000  FT 

-0.382974 

0 

0 

KDOV  SKY  (OBSCURED) 

0.1530178 

1 

0.1530178 

KDOV  VIS  <  4800  M 

-0.139053 

0 

0 

Forecast  Probability: 

0.7690357 

place  on  24  Jan  01  at  0627Z  at  Dover  AEB  and  the  forecast  is  for  ceiling  below  200  ft. 
Eigure  13  illustrates  the  two  probability  forecasts  for  the  six  valid  times  along  with  the 
actual  observed  event.  Notice  at  the  one-hour  point,  the  observations-based  network 
produces  a  forecast  of  over  80%  probability  of  occurrence,  while  conditional  climatology 
has  probabilities  in  the  60-65%  range.  As  the  low  ceiling  lifts,  slightly  after  0800Z,  the 
two-hour  observations-based  forecast  accurately  detects  dissipation,  generating  a 
probability  forecast  of  17%  compared  to  33%  for  conditional  climatology.  The 
observations  forecast  remains  in  sync  with  the  actual  observation  at  the  three-hour  point. 
The  four-hour  point  shows  both  methods  producing  forecast  probabilities  of  25%,  still 
within  the  good  forecast  range.  Einally,  the  observations-based  network  accurately 
forecasts  the  last  two  hours,  while  conditional  climatology  lags  behind  in  its  prediction. 

The  second  example  (Eig.  14)  is  from  26  Oct  00,  0500Z  at  McGuire  AEB,  with 
visibility  less  than  0.5  mi  as  the  forecasted  predictand.  This  situation  persists  for  over  six 
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KDOV  Ceiling  Less  than  200  ft  Examjde 


-OBS 

-CC 

-A- 

-  Actual 

FIG.  13.  Forecast  Comparison  for  Dover  AFB.  Each  forecast  method 
is  plotted  for  a  six  hour  forecast  beginning  at  0627Z,  24  Jan  01  at  Dover 
AFB.  Note  that  the  observations-based  network  predicts  dissipation 
faster  and  more  accurately  than  conditional  climatology. 


hours  and  is  forecasted  more  accurately  by  the  observations-based  network.  The  one- 
hour  forecast  is  accurate  for  both  forecast  techniques,  with  a  92%  probability  for  the 
observations-based  network  and  an  81%  probability  for  conditional  climatology.  The 
second  hour  illustrates  the  strength  of  the  observations-based  forecast,  where  conditional 
climatology  almost  always  decreases  probabilities  with  time;  the  observations  network 
adjusts  and  produces  a  more  accurate  forecast  based  solely  on  the  current  conditions. 

Thus  the  two-hour  forecast  is  actually  more  accurate  than  the  one-hour  forecast  and 
significantly  better  than  conditional  climatology.  Although  the  probabilities  decrease  as 
time  goes  on,  the  observations-based  method  continues  to  produce  a  superior  forecast. 

As  far  as  five  hours  out,  this  method  is  still  forecasting  a  64%  probability  of  occurrence, 
which  may  be  below  the  good  forecast  criteria,  but  is  still  a  solid  forecast  for  five  hours  in 
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the  future  based  solely  on  the  current  conditions  and  is  significantly  better  than 
conditional  climatology. 

Overall  these  two  examples  illustrate  the  effectiveness  of  an  observations-based 
system  to  accurately  forecast  timing  of  dissipation,  as  well  as  event  duration. 
Conclusions  arising  from  this  research  along  with  recommendations  on  implementation 
and  future  work  follow. 


KWRI  Visibility  Less  than  0,5  mi  Example 


-•-OBS 
«-CC 
.A— Actual 


FIG.  14.  Forecast  Comparison  for  McGuire  AFB.  Each  forecast 
method  is  plotted  for  a  six  hour  forecast  beginning  at  0500Z,  26  Oct  00 
at  McGuire  AFB.  Note  that  the  observations-based  network  more 
accurately  predicts  the  continued  presence  of  the  low  ceiling  through  the 
duration  of  the  event. 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions 

The  goal  of  this  research  is  to  develop  a  forecast  tool  for  the  15  OWS  to  predict 
dissipation  of  fog  more  accurately  than  their  current  methodology,  conditional 
climatology.  To  accomplish  this,  36  probabilistic  forecast  equations  are  developed,  12 
for  each  of  the  three  main  airfields  along  the  east  coast,  Dover  AFB,  McGuire  AFB,  and 
Andrews  AFB.  Short-term  forecast  equations  are  developed  at  hourly  intervals  from  0-6 
hours  beyond  the  initial  time  for  each  of  two  predictors,  ceiling  less  than  200  ft  and 
visibility  less  than  0.5  mi. 

The  overall  results  of  this  research  are  mixed.  Of  the  36  predictive  equations,  16 
perform  better  statistically  than  conditional  climatology  when  compared  against  an 
independent  sample  of  data.  Of  the  forecasts  that  are  better,  a  17.6%  average 
improvement  is  made  over  conditional  climatology.  Some  of  the  improved  forecasts  are 
the  ceiling  less  than  200  ft  criteria  at  both  Dover  AFB  and  McGuire  AFB  where  all  6 
hourly  forecasts  had  positive  skill  scores.  Also  notable  are  the  McGuire  AFB  visibility 
forecasts,  which  show  a  positive  average  skill  score.  The  two  most  successful  forecasts, 
for  ceilings  less  than  200  ft  at  Dover  AFB  and  McGuire  AFB  both  show  large  increases 
in  skill  score,  a  high  number  of  good  forecasts,  and  few  bad  forecasts. 
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In  the  cases  where  the  observations-based  network  is  superior,  the  predictive 
equations  have  similar  characteristics.  First,  is  geographic  influence.  The  early  forecasts 
have  little  dependence  on  the  geographic  location  of  the  network  station  selected.  As  the 
forecast  lead  time  increases  dependence  on  predictors  from  the  southwest  especially 
dominate  the  equations.  Second,  many  of  the  successful  longer  term  forecasts  (5  or  6- 
hours)  contain  climatological  parameters.  A  strong  dependence  on  wind  predictors,  both 
at  the  surface  and  upper  levels  is  evident  throughout  the  equations. 

With  these  similarities  examined,  it  is  noted  that  each  of  these  forecasts  is 
independently  produced  and  validated  for  each  station  at  each  hour.  Although  there 
appear  to  be  trends,  switching  a  predictor  from  one  model  to  another  will  not  produce  a 
more  accurate  forecast.  Just  because  one  predictor  works  at  McGuire  AFB  does  not 
mean  it  will  be  effective  at  Dover  AFB. 

The  most  disappointing  conclusions  of  this  study  are  the  Andrews  AFB  forecasts, 
which  show  large  negative  skill  scores.  A  majority  of  the  forecast  times  for  Andrews 
AFB  fail  the  tests  of  normality  and/or  heteroscidacity  required  to  meet  the  assumptions 
necessary  to  employ  linear  regression,  indicating  that  a  linear  model  is  not  the  best 
representation  of  this  data  set.  Problems  contributing  to  this  include  the  Andrews  AFB 
data  set  being  smaller  than  the  others,  especially  in  the  veriflcation  data  set  where 
missing  observations  dominate.  As  a  result,  the  Andrews  AFB  data  set  is  unable  to 
account  for  variability  in  surface  observations  (in  fog  cases),  which  prevents  the 
production  of  an  accurate  forecast.  More  research,  with  a  larger,  more  complete  data  set 
may  correct  these  forecasts  for  operational  implementation. 
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Despite  forecast  scores  below  conditional  climatology,  the  negative  forecasts  do 
show  some  advantages.  First,  the  number  of  forecasts  with  absolute  error  less  than  30 
percent  (e.g.  good  forecasts)  often  exceeds  the  number  of  good  forecasts  produced  by 
conditional  climatology.  Second,  the  MSEs  in  these  cases  are  typically,  with  a  few 
exceptions,  low  and  very  competitive  with  the  climatology  forecasts.  Many  of  the  errors 
and  decreased  forecasting  skill  from  the  poor  forecasts  in  the  observations-based  network 
are  similar  and  shown  to  have  bias.  The  most  common  bias  among  the  forecasts  is  an 
overestimate  of  event  occurrence  when  the  fog  has  already  dissipated.  In  the  later 
forecast  hours  of  the  forecast  this  is  especially  significant  and  results  in  the  higher 
degradation  of  skill  scores. 

A  final  success  noted  is  the  overall  application  of  the  methodology  to  new 
“extreme”  predictand  values,  representing  the  worst  cases  of  fog.  The  lowest  values 
studied  before  now  were  500  ft  ceilings  and  1  mi  visibility  (Vislocky  and  Fritsch  1997). 
This  is  encouraging  as  this  allows  the  method  to  have  significant  impact  on  operations  in 
the  most  marginal  weather  conditions.  Recommendations  for  implementation  of  this 
technique  by  the  15  OWS  and  possible  future  research  to  improve  the  shortfalls  follow. 

5.2  Recommendations  for  the  15  OWS 

It  is  recommended  that  the  15  OWS  implement  this  program  for  the  following 
cases:  KDOV  ceiling  less  than  200  ft,  KWRI  ceiling  less  than  200  ft,  and  KWRI  visibility 
less  than  0.5  mi.  For  both  of  the  ceiling  forecasts,  all  six  hourly  forecasts  show 
significant  improvement  over  conditional  climatology  and  are  the  superior  forecast 
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model.  The  six  hourly  visibility  foreeasts  for  KWRI  only  have  three  hours  that  are 
eonsistently  better  than  eonditional  elimatology;  however,  it  has  two  advantages.  First, 
the  earlier  hours  eonsistently  produee  the  better  foreeasts;  only  the  first  hour  had  a 
negative  skill  seore.  Seeond,  the  first  hour  had  a  distinet  advantage  in  produeing  good 
foreeasts— foreeasts  with  an  absolute  error  of  less  than  30%.  Therefore,  this  model  ean  be 
used  sueeessfully  in  an  operational  environment. 

Implementation  of  the  other  three  models  eould  be  used  on  a  trial  basis.  Although 
the  verification  results  are  not  optimal,  the  technique  has  some  merit  in  these  cases  based 
on  both  previous  work  and  this  current  research.  Implementation  of  these  methods  on  a 
trial  basis  could  produce  better  results  than  in  the  confined  scope  of  this  research. 
Realizing  the  shortfalls  of  each  particular  model  is  crucial  for  successful  implementation. 
For  example,  the  Dover  AFB  visibility  forecasts,  which  all  have  negative  skill  scores  do  a 
solid  job  of  forecasting  occurrences  in  the  short-term,  however,  the  oveforecasting 
tendency  is  tied  to  the  reliance  on  southwesterly  winds  to  advect  the  fog  from  the  area. 
Analysis  of  the  presence  or  absence  of  these  conditions  could  lead  the  forecaster  to  an 
accurate  dissipation  forecast.  Using  these  models  on  any  type  of  basis  could  provide 
valuable  information  on  the  ability  of  the  model  to  produce  successful  forecasts. 

5.3  Recommendations  for  Future  Research 

The  observations-based  method  is  successful  for  many  areas,  including  the  area 
currently  researched  for  the  15  OWS.  This  researcher  believes  that  with  some  tweaking 
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of  the  system,  this  methodology  eould  be  applied  across  the  AF.  There  are  factors  that 
must  be  accounted  for  in  future  research  to  make  this  work  operationally. 

First,  forecasts  should  be  broken  up  temporally.  Forecast  increments  of  2-4  hours 
should  be  looked  at  individually  in  order  to  focus  on  the  exact  physical  mechanisms  that 
dissipate  the  fog.  This  allows  diurnal  cycles  to  be  taken  into  account,  capitalizing  on  the 
advantages  conditional  climatology  offers.  Distinct  time  periods  would  be  more 
representative  of  different  observed  parameters  predicting  dissipation. 

Vislocky  and  Fritsch  (1997)  suggest  inclusion  of  satellite  imagery  into  the  system. 
This  is  a  key  element  to  determine  spatially  how  large  the  extent  of  the  fog  area  is. 
Although  the  statistical  method  is  cut  and  dry,  the  input  of  both  the  spatial  extent  of  the 
fog  coverage,  as  well  as  how  close  the  airfield  of  interest  is  to  the  outer  edges  of  the  fog 
(as  suggested  by  Gurka  1978)  would  add  considerable  predictive  ability  to  the  model. 

Additionally,  further  research  is  necessary  to  apply  this  technique  throughout  the 
remainder  of  the  15  OWS  AOR.  The  original  intent  of  the  research  was  to  apply  this 
methodology  to  all  1 1  airfields  throughout  the  15  OWS.  With  some  modification  to  the 
methodology,  this  could  be  a  successful  application  for  any  site  within  the  AOR.  It  is 
envisioned  that  this  methodology,  with  a  significant  amount  of  additional  work,  could  be 
applied  to  airfields  throughout  the  AF  as  an  effective  alternative  to  conditional 
climatology. 

This  study’s  observations-based  methodology  could  benefit  significantly  from  the 
inclusion  of  mesonet  data  in  regions  where  available.  States  such  as  Oklahoma  and 
Colorado  currently  have  mesonets.  This  spatially  denser  observing  platform  could 
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provide  more  insight,  data,  and  information  on  the  dissipation  process  and  as  a  result 
produce  more  accurate  forecasts. 

Lastly,  another  process  that  could  produce  accurate  forecast  equations  is  the  use 
of  a  logistic  regression  model  as  opposed  to  the  linear  model.  The  advantage  of  the 
logistic  model  is  that  it  always  produces  a  value  between  0  and  1  and  is  able  to  reduce 
forecast  error  (Hilltker  and  Fritsch  1999).  A  comparison  between  linear  and  logistical 
regressions  could  demonstrate  the  advantages  of  the  logistical  regression  and  further 
improve  forecasts. 
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Appendix:  Surface  Observing  Networks 


This  appendix  outlines  the  40  surface  observing  stations  used  to  produce  the 
dissipation  forecasts  for  the  three  airfields.  The  10,  19,  25,  30,  35,  and  40  station 
networks  used  to  produce  the  1-6  hour  forecasts  are  listed. 
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Dover  AFB 


TABLE  A1 .  Surface  Network  Observing  Stations  for  Dover  AFB. 


Observing  Station 

ICAO 

EAT 

LON 

Elevation  (m) 

10  Station  Network 

Dover  AFB  (DE) 

KDOV 

39  08 

75  28 

9 

Georgetown  Suffox  County  (DE) 

KGED 

38  41  24 

75  21  45 

15 

Wilmington  (DE) 

KILG 

39  40  22 

75  36  03 

24 

Wildwood  (NJ) 

KWWD 

39  01 

74  55 

7 

Millville  (NJ) 

KMIV 

39  21  58 

75  04  42 

26 

Philadelphia  (PA) 

KPHL 

39  52  06 

75  14  37 

6 

Philadelphia  NE  Airport  (PA) 

KPNE 

40  04  44 

75  00  49 

36 

Salsbury  (MD) 

KSBY 

38  20  21 

75  30  15 

15 

Atlantic  City  (NJ) 

KACY 

39  27  53 

074  35  12 

23 

Baltimore,  Washington  International  (MD) 

KBWI 

39  10  00 

076  41  00 

44 

19  Station  Network 

Baltimore,  Martin  (MD) 

KMTN 

39  20 

076  25 

7 

Mount  Holley  (NJ) 

KVAY 

39  56  26 

074  50  28 

16 

Patuxent  River  NAS  (MD) 

KNHK 

38  16  43 

076  24  50 

12 

Andrews  AFB  (MD) 

KADW 

38  49 

076  51 

86 

Wallops  Island  (VA) 

KWAL 

37  56  26 

075  27  47 

12 

McGuire  AFB  (NJ) 

KWRI 

40  01 

074  36 

41 

Trenton,  Mercer  County  Airport  (NJ) 

KTTN 

40  16  35 

074  48  59 

64 

Reading  Regional  Airport  (PA) 

KRDG 

40  22  24 

075  57  34 

105 

Lancaster  Airport  (PA) 

KENS 

40  07  13 

076  17  40 

123 

25  Station  Network 

Allentown  (PA) 

KABE 

40  39  03 

075  26  57 

120 

Fort  Belvoir  (VA) 

KDAA 

38  43 

077  11 

21 

Washington,  Reagan  National  (VA) 

KDCA 

38  50  54 

077  02  03 

4 

Harrisburg,  Capitol  City  Airport  (PA) 

KCXY 

40  13  13 

076  51  14 

105 

Harrisburg  International  Airport  (PA) 

KMDT 

40  1 1  46 

076  46  23 

94 

York  Airport  (PA) 

KTHV 

39  55  22 

076  52  41 

146 

30  Station  Network 

Belmar-Farmsdale  (NJ) 

KBLM 

40  11 

074  08 

48 

Melfa/ Accomack  Airport  (VA) 

KMFV 

37  39 

075  46 

15 

Quantico,  MCAF  (VA) 

KNYG 

38  30  45 

077  17  30 

3 

Washington-Dullus  (VA) 

KIAD 

38  56  05 

077  26  51 

95 

Shannon  Airport  (VA) 

KEZF 

38  16 

077  27 

26 

35  Station  Network 

Culpeper  County  Airport  (VA) 

KCJR 

38  31  36 

077  51  32 

95 

Newerk  International  Airport  (NJ) 

KEWR 

40  40  57 

074  10  10 

5 

Caldwell,  Exxex  Coimty  Airport  (NJ) 

KCDW 

40  52  35 

074  16  59 

52 

Charlottesville-Abermarle  Airport  (VA) 

KCHO 

30  08  18 

078  27  21 

195 

Richmond  International  Airport  (VA) 

KRIC 

37  30  40 

077  19  24 

51 

40  Station  Network 

Farmingdale  Republic  Airport  (NY) 

KFRG 

40  44  03 

073  25  01 

24 

New  York  City,  JFK  (NY) 

KJFK 

40  38  19 

073  45  44 

3 

Langly  AFB  (VA) 

KLFI 

37  05 

076  21 

6 

New  York,  LaGuardia  (NY) 

KLGA 

40  46  45 

073  52  48 

6 

Petersburg  (VA) 

KPTB 

37  11 

077  31 

59 
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McGuire  AFB 


TABLE  A2.  Surface  Network  Observing  Stations  for  McGuire  AFB. 


Observing  Station 

ICAO 

LAT 

LON 

Elevation  (m) 

10  Station  Network 

McGuire  AFB  (NJ) 

KWRI 

40  01 

074  36 

41 

Mount  Holley  (NJ) 

KVAY 

39  56  26 

074  50  28 

16 

Belmar-Farmsdale  (NJ) 

KBLM 

40  11 

074  08 

48 

Trenton,  Mercer  County  Airport  (NJ) 

KTTN 

40  16  35 

074  48  59 

64 

Philadelphia  International  Airport  (PA) 

KPHL 

39  52  06 

075  13  52 

6 

Philadelphia,  NE  Philadelphia  Airport  (PA) 

KPNE 

40  04  44 

075  00  49 

36 

Atlantic  City  (NJ) 

KACY 

39  27  53 

074  35  12 

23 

Millville  Municipal  Airport  (NJ) 

KMIV 

39  21  58 

075  04  42 

26 

Wilmington  (DE) 

KILG 

39  40  22 

075  36  03 

24 

Allentown  (PA) 

KABE 

40  39  03 

075  26  57 

120 

19  Station  Network 

Newerk  International  Airport  (NJ) 

KEWR 

40  40  57 

074  10  10 

5 

Caldwell,  Exxex  County  Airport  (NJ) 

KCDW 

40  52  35 

074  16  59 

52 

New  York  City,  JFK  (NY) 

KJFK 

40  38  19 

073  45  44 

3 

New  York,  LaGuardia  (NY) 

KLGA 

40  46  45 

073  52  48 

6 

Wildwood  ASOS  (NJ) 

KWWD 

39  01 

074  55 

2 

Teterboro  Airport  (NJ) 

KTEB 

40  51  32 

074  03  24 

2 

Islip,  Long  Island  (NY) 

KISP 

40  47  38 

073  06  06 

30 

Farmingdale  Republic  Airport  (NY) 

KFRG 

40  44  03 

073  25  01 

24 

Dover  AFB  (DE) 

KDOV 

39  08 

075  28 

9 

25  Station  Network 

Reading  Regional  Airport  (PA) 

KRDG 

40  22  24 

075  57  34 

105 

Lancaster  Airport  (PA) 

KLNS 

40  07  13 

076  17  40 

123 

Georgetown,  Suxxex  County  Airport  (DE) 

KGED 

38  41  24 

075  21  45 

15 

Montgomery,  Orange  County  Airport  (NY) 

KMGJ 

41  30  33 

074  15  54 

111 

Newburgh/Stewart  (NY) 

KSWF 

41  30 

074  06 

150 

Wilkes-Barre-Scranton  International  (PA) 

KAVP 

41  20  20 

075  43  36 

293 

30  Station  Network 

Bridgeport  (CT) 

KBDR 

41  09  30 

073  07  44 

3 

Danbury  Municipal  Airport  (CT) 

KDXR 

41  22  18 

073  29  04 

139 

Westhampton  Beach  (NY) 

KFOK 

40  51  03 

072  37  14 

20 

Baltimore,  Martin  (MD) 

KMTN 

39  20 

076  25 

7 

Baltimore,  Washington  International  (MD) 

KBWI 

39  10  00 

076  41  00 

44 

35  Station  Network 

Monticello  (NY) 

KMSV 

41  42 

074  48 

428 

Poughkeepsie  (NY) 

KPOU 

41  37  32 

073  52  55 

50 

Salsbury  (MD) 

KSBY 

38  20  21 

75  30  15 

15 

Andrews  AFB  (MD) 

KADW 

38  49 

076  51 

86 

Washington,  Reagan  National  (VA) 

KDCA 

38  50  54 

077  02  03 

4 

40  Station  Network 

Fort  Belvoir  (VA) 

KDAA 

38  43 

077  11 

21 

Washington-Dullus  (VA) 

KIAD 

38  56  05 

077  26  51 

95 

Harrisbing,  Capitol  City  Airport  (PA) 

KCXY 

40  13  13 

076  51  14 

105 

Harrisbing  International  Airport  (PA) 

KMDT 

40  1 1  46 

076  46  23 

94 

York  Airport  (PA) 

KTHV 

39  55  22 

076  52  41 

146 

74 


Andrews  AFB 


TABLE  A3.  Surface  Network  Observing  Stations  for  Andrews  AFB. 


Observing  Station 

ICAO 

EAT 

LON 

Elevation  (m) 

10  Station  Network 

Andrews  AFB  (MD) 

KADW 

38  49 

076  51 

86 

Washington,  Reagan  National  (VA) 

KDCA 

38  50  54 

077  02  03 

4 

Fort  Belvoir  (VA) 

KDAA 

38  43 

077  11 

21 

Manassas  Municipal  (VA) 

KHEF 

38  43 

077  31 

59 

Quantico  MCAF  (VA) 

KNYG 

38  30  45 

077  17  30 

3 

Washington-Dullus  (VA) 

KIAD 

38  56  05 

077  26  51 

95 

Baltimore,  Washington  International  (MD) 

KBWI 

39  10  00 

076  41  00 

44 

Baltimore,  Martin  (MD) 

KMTN 

39  20 

076  25 

7 

Shannon  Airport  (VA) 

KEZF 

38  16 

077  27 

26 

Hagerstown  (MD) 

KHGR 

39  42  21 

077  43  48 

214 

19  Station  Network 

Culpeper  County  Airport  (VA) 

KCJR 

38  31  36 

077  51  32 

95 

Winchester  Regional  (VA) 

KOKV 

39  09 

078  09 

222 

Martinsburg  (WV) 

KMRB 

39  24  14 

077  58  30 

169 

York  Airport  (PA) 

KTHV 

39  55  22 

076  52  41 

146 

Patuxent  River  NAS  (MD) 

KNHK 

38  16  43 

076  24  50 

12 

Richmond  International  Airport  (VA) 

KRIC 

37  30  40 

077  19  24 

51 

Dover  AFB  (DE) 

KDOV 

39  08 

75  28 

9 

Lancaster  Airport  (PA) 

KENS 

40  07  13 

076  17  40 

123 

Harrisburg,  Capitol  City  Airport  (PA) 

KCXY 

40  13  13 

076  51  14 

105 

25  Station  Network 

Harrisburg  International  Airport  (PA) 

KMDT 

40  1 1  46 

076  46  23 

94 

Georgetown  Suffox  County  (DE) 

KGED 

38  41  24 

75  21  45 

15 

Wilmington  (DE) 

KILG 

39  40  22 

75  36  03 

24 

Salsbury  (MD) 

KSBY 

38  20  21 

75  30  15 

15 

Wallops  Island  (VA) 

KWAL 

37  56  26 

075  27  47 

12 

Melfa/ Accomack  Airport  (VA) 

KMFV 

37  39 

075  46 

15 

30  Station  Network 

Petersburg  (VA) 

KPTB 

37  11 

077  31 

59 

Charlottesville-Abermarle  Airport  (VA) 

KCHO 

30  08  18 

078  27  21 

195 

Reading  Regional  Airport  (PA) 

KRDG 

40  22  24 

075  57  34 

105 

Philadelphia  International  Airport  (PA) 

KPHL 

39  52  06 

075  13  52 

6 

Philadelphia  NE  Airport  (PA) 

KPNE 

40  04  44 

75  00  49 

36 

35  Station  Network 

Wildwood  ASOS  (NJ) 

KWWD 

39  01 

074  55 

2 

Millville  Municipal  Airport  (NJ) 

KMIV 

39  21  58 

075  04  42 

26 

Langly  AFB  (VA) 

KLFI 

37  05 

076  21 

6 

Staunton/Shenandoah  (VA) 

KSHD 

38  16 

078  54 

366 

Atlantic  City  (NJ) 

KACY 

39  27  53 

74  34  02 

23 

40  Station  Network 

Trenton,  Mercer  County  Airport  (NJ) 

KTTN 

40  16  35 

074  48  59 

64 

Newerk  International  Airport  (NJ) 

KEWR 

40  40  57 

074  10  10 

5 

Belmar-Farmsdale  (NJ) 

KBLM 

40  11 

074  08 

48 

Allentown  (PA) 

KABE 

40  39  03 

075  26  57 

120 

McGuire  AFB  (NJ) 

KWRI 

40  01 

074  36 

41 

75 
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